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Summary— In a previous paper [5] the authors have discussed the development of a correlation 
for the prediction of heat transfer coefficients of pseudo-plastic fluids. For the application 
of the correlation to any particular problem, it is necessary to have viscosity data over the 
range of conditions involved. The present paper discusses the method employed for obtaining 
the viscosity data used in the correlation of heat transfer coefficients. 

Three plastic (pseudo-) fluids were investigated covering a range of viscosities from one to 
thirteen centipoises. The pseudo-plastic liquids used were the aqueous solution of sodium carboxy- 
methyl cellulose, polyvinyl aleohol, and sodium polymethacrylate. The modified Ostwald 
viscometer was designed and used in this work. Previous workers have used either pipe visco- 
meters or the rotating cylinder type, both of which are cumbersome and expensive to install. 
The modified Ostwald viscometer maintains the condition of tubular flow and conducts the 
measurement at very low rates of shear. Another advantage of the modified Ostwald viscometer 
is that it enables viscosity measurements to be made at point values of the pressure drop across 
the capillary tube. This is important for pseudo-plastic fluids when the viscosity is dependent 
upon the rate of flow. 

The results of viscosity measurement can be directly applicable to the design of heat 
exchangers and pressure drop calculations for the specification of the pump handling the pseudo- 
plastic fluids. It is probable that a similar treatment can be applied to those fluids exhibiting 


general non-Newtonian behaviour. 


Résumé— Dans un mémoire antérieur, les auteurs ont discuté létablissement d'une relation 
permettant la prévision des coefficients de transmission thermique des fluides pseudo-plastiques. 
Pour appliquer cette relation & un probleme particulier, il faut connaitre les valeurs de la viscosité 
dans l’étendue du domaine d'utilisation considéré. Les auteurs discutent ici la méthode qu’‘ils 
emploient pour obtenir les valeurs de la viscosité utilisées dans la corrélation des coefficients 
de transmission thermique. 

Ils ont étudié trois fluides pseudo-plastiques couvrant des valeurs de la viscosité comprises 
entre une et treize centipoises. Les liquides pseudo-plastiques étaient une solution aqueuse 
d’oxyméthyl-cellulose-carbonate de sodium, d’alcool polyvinylique et de polyméthacrylate de 
sodium. Ils ont utilisé un viscosimétre d’'OstTwaLp convenablement modifié, tandis que les auteurs 
antérieurs avaient employé des \ iscosimétres a tube ou a cylindre rotatif, appareils encombrants et 
d'installation coiteuse. Le viscosimétre d’'OstwaLp modifié assure un état d’écoulement laminaire 
et permet la mesure avec des tensions de cisaillement trés faibles. L’appareil présente l'avantage 
supplémentaire de permettre des mesures de viscosité pour des valeurs parfaitement définies 
de la perte de charge a travers le capillaire : ceci est important pour les liquides pseudo-plastiques 
dont la viscosité varie avec la vitesse d’écoulement. 

Les résultats de ces mesures sont directement applicables au dessin des échangeurs de 
chaleur et aux calculs des pertes de charges indispensables pour le choix des pompes pour la 
manutention des liquides pseudo-plastiques. Il parait vraisemblable que le méme mode de traite- 
ment serait applicable & d'autres fluides non-Newtoniens. 
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INTRODUCTION 


In a previous paper (5) the authors have developed 
a correlation for the prediction of heat transfer 
coefficients of pseudo-plastic fluids. The concept 
of heat transfer coefficient the 


is based upon 


assumption of a lamina flow adjacent to the 


wall. For both Newtonian and Non-Newtonian 
fluids, a knowledge of viscosity characteristics 
of the fluid is a prerequisite to the development 
of a correlation of heat transfer coefficient data. 
The viscosity of a fluid is defined as the rate 
of a shear stress (Sg) to the rate of shearing (R,) 
perpendicular to the direction of the shear stress. 
The rate of shearing, R,, is expressed in terms 
of velocity gradient. Mathematically, the viscosity 
is given by : 
p = Sg/R, (1) 
This ratio is a constant independent of shear 
rate for a majority of fluids which are known as 
Fluids. At 


the shear rate is equal to zero. For a flow through 


Newtonian zero stress of shear, 


a channe] of circular cross section, the viscosity 
of a fluid at the wall is given by 


PRg/2N (2 
rk? 
where 


PRg/QN = S (3) 


= R, (4) 
the 
wall can be evaluated from the volumetric rate 
of flow of fluids and the diameter of the circular 
The the 
evaluated from the length and diameter of the 


Eq. (2) indicates that the shear rate at 


channel shearing stress at wall is 
channel with the data of pressure drop through 
the of flow. the 
viscosity fluid is a constant 
for a given fluid at a given temperature and 
independent of shear stress, the visocity at the 


wall is the same as the viscosity of the fluid at 


channel in a direction Since 


of a Newtonian 


any point in an isothermal flow. 

On the other hand, if the ratio of shear stress to 
shear rate is not a constant, the fluid is called non- 
Newtonian. Pseudo-plastic fluids come into this 
general catagory. The viscosity of a pseudo- 


plastic fluid has a finite alue at zero rate of shear, 
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falling off asymptotically to a lower fixed value 
as the rate of shear increases to infinity. There 
is no time dependence of the viscosity on the 
condition of shear. Examples of these fluids are 
many suspensions, latices, emulsions and solutions 
of larger molecules. 

The second class of fluids deviating from 
behaviour are those known as 
These 


evidence of flow until a certain minimum value 


Newtonian 


“ Bingham bodies.” materials show no 


of shear stress is achieved. Beyond this point 


commence to flow and exhibit simple 


they 
Newtonian behaviour providing the shear stress 
is maintained in excess of the yield point value. 
Such a 
plastic. In 


material is sometimes termed an ideal 


practice many materials behave 
as Bingham bodies insofar as they show yield 
points, but thereafter they behave as pseudo- 
plastic fluids rather than as simple Newtonian 
fluids. These fluids are termed by Hepstrom [8] 
to be general non-Newtonian fluids. Fig. 1 shows 
graphically the shearing stress vs. rate of shear 
relation for the various types of fluid discussed. 

The 


pseudo-plastic fluids, but it ts probable that a 


present work is confined primarily to 


RATE OF SHEAR 


A = SIMPLE NEWTONIAN FLUID 


= PSEUDO PLASTIC FLUID 


= BINCHAM BODY 


D= CENERAL NON NEWTONIAN FLUIO 


Fig. 1. Shear rate and stress diagram for fluids :— 
A = Simple Newtonian fluid. 

B = Pseudo-plastic fluid 

Cc Binham body 


D General non-Newtonian fluid. 
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similar treatment can be applied to those fluids 
exhibiting general non-Newtonian behaviour. 
There are several terms with regard to the 
viscometric behaviour of a pseudo-plastic fluid. 
For the isothermal flow of pseudo-plastic fluid, 
Wiwnpinc has successfully used in his correlation 
[16-18] “ apparent viscosity’ which was defined 
by (2) (12] (16) [18] 
= Sg/R, (5) 


The rate of shear, R,, and shearing stress Sg 
are calculated the same was as for a Newtonian 
fluid by Eqs. (3) and (4). 
shear computed from the volume by Eq. (4) 


However, the rate of 


can no longer be the true rate of shear at the 
wall because of its dependence upon the shear 
stress. Another term is called the differential 
viscosity, which is defined as the slope at any 
point of the curve of shear versus shear rate. 
Mathematically, it is given by 
dSg 


1 (6) 
aR, 


Fig. 2 shows the shearing stress vs. shearing 
rate curve of a_ typical pseudo-plastic fluid. 
It will be seen that the differential 


at any point X on the curve is the slope at point 


viscosity 


X, while the apparent viscosity is the ratio ii 


R, 
B= CURVE FOR PSEUDO PLASTIC FLUID 


De CURVE FOR CENERAL NON NEWTONIAN FLUID 


Fig. 2. Sheer rate and stress diagram for non-Newtonian 
fluids. 
B = Curve for pseudo-plastic fluid. 
D = Curve for general non-Newtonian fluid. 
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Although the apparent and differential viscosi- 


ties are empirical quantities, they represent 
unique measures of the viscometric behaviour 
of a pseudo-plastic fluid at any selected condition 
of shear, and as such they are satisfactory for 
most engineering purposes. 

There have been many attempts [1, 4, 11, 12, 
13, 14, 15, 16] to derive an expression in relating 
the shear stress to the rate of shear for pseudo- 
plastic fluids. These expressions are known as 
14], the coefficient 


of structural stability method [14], the direct 


the power series method (13, 


integration of tubular flow equation [11], and 
WILLIAMSONS equation [15]. Using a hyperbolic 
function of shear stress, Brown and BuRRIDGE 
[4] developed a relation which was too cumber- 
practical application. 


some to permit § any 


equation [15], which the 
simplest relation between the shear stress and 
shear rate of a pseudo-plastic fluid, was used by 
Winpinc et al. in the treatment of their data 
[16, 17, 18]. The equation is of the following 


form : 


(7) 


The constants in the equation, S, g, a, uw. are 
functions of the system and temperature and 
not of rate of shear nor the dimensions of the 
tube. 

As a consequence of the various methods of 
defining the coefficient of pseudo-plastic visocosity, 
it follows that the Reynolds number, depending 
as it does on viscosity, may also be expressed in 
several ways for a pseudo-plastic fluid. WInDING 
[16] recommended the use of the apparent 


viscosity, giving an apparent Reynolds number, 


(8) 


While the apparent viscosity is not a constant 
quantity across the section of a pseudo-plastic 
fluid in tubular flow, its use may be justified 
on the same basis as the use of the mass velocity, 
which also quantity, 


represents an average 


in the calculation of the Reynolds number 
[16]. In deriving the pseudo-plastic analogue 


of the Fanning friction factor, Winpino et al. 
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(18) found that data on GR-S latices were 
correlated better by the use of apparent viscosity 
at zero rate of shear. 

In work on slurries behaving essentially as 
“ Bingham bodies,”” [3] and Hepstrom 
[8] derived the apparent Reynolds number in 
terms of effective viscosity of slurry. The effective 
viscosity is mathematically analogous to the 
apparent viscosity at of shear, 
#..» of a pseudo-plastic fluid. In fact, Winpinc 
[17, 18] used the apparent viscosity at infinite 
rate of shear in the correlation of some data 
particularly on heat transfer. 


infinite rate 


TO AC MAINS 


Fig. 3. General arrangement drawing of the viscosity 
apparatus. 


EXPERIMENTAL MetTHop 
A. Description of the apparatus. 
A schematic drawing of the general arrangement 


of the viscosity apparatus is presented in Fig. 3. 
Detailed dimensional drawings of the modified 


Ostwald viscometer used are shown in Fig. 4. 
A glass thermostat bath (A) was filled with water 
and temperature controlled to + 0-05°C, the 
temperature of the bath being indicated by the 
precision grade mercury in glass thermometer 
(F) graduated in 0-1°C, The marine type impeller 
(E) assured vigorous agitation of the water in the 
bath. 

In the thermostat bath were immersed the 
viscometer (G) and the air bubbler (H). The 
water level in the bath was maintained just 
below the cocks (K,, K,, and K,) on the viscometer. 
The exit tubes (J, and J,) of the viscometer 


35, 80 (40 


“| CAP 


LLARY. 


ALL DIMENSIONS IN CMS 


Fig. 4. Modified Ostwald viscometer. 


were connected via pressure rubber tubing 
to the air feed and exit manifold system com- 
prising cocks K,, K,, K,, K,, and L. This 
manifold system provided the flexibility of air 
supply, and of venting to atmospheric pressure, 
which was necessary for the development of a 
convenient experimental procedure. This mani- 
fold also had a connection via cock K, to the 
mercury manometer (M), or via cock K, to the 
water manometer (N). 

Filtered compressed air was supplied from 
an 80 p.s.i.g. main via a reducing value (P). 
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Suitable adjustment of this reducing valve 
enabled a constant pressure on the exit side of 
the valve to be maintained at any value greater 
than 3 cm Hg gauge. For lower exit pressures the 
operation of the reducing valve became erratic. 
Constant pressures below 83cm Hg gauge were 
obtained by adjusting the valve to give an exit 
pressure of 8cmHg gauge and introducing 
a leak to atmosphere on the exit side of the valve 
by adjusting the screw clip (L). The bubbler (H) 
was necessary to saturate the feed air at the 
operating temperature of thermostat in order 
to avoid evaporation of the sample in the visco- 
meter. The exit pipe from the bubbler to the 
and K, was lagged. 


T-joint leading to cocks K, 


B. Experimental procedure 


The “ Ostwald viscometer,” used in this work 
differed from those used ordinarily in that the 
bulbs were designed to hold a sample of approxi- 
mately 100 ml. This large sample size minimized 
the effect of any errors in measuring the time of 
passage of the sample when flowing through the 
capillary at high rate. Such measurements would 
not be necessary in the case of a simple Newtonian 
fluid but becomes essential with a pseudo-plastic 
fluid if its viscometric behaviour is to be 
established over a reasonably wide range of 
rates of shear. Four such viscometers were 
constructed. These viscometers were of essentially 
identical external dimensions and differed only 
in the bore of the capillary tube. The bore 
ranged from about 0-075cm to 0-2cm. The 
details of the viscometer dimensions are given 
in Table 1. 

The capillaries were specially selected to be 
of uniform bore to within + 0-4% throughout 
their length, and the sample volumes were 
measured to + 0-1 ml. 

It was necessary to obtain data from visco- 
meters over a range of capillary diameters in 
order to determine the significance of the wall 
slippage effect in relation to capillary diameter 
[16]. The viscometer was filled with approxi- 
mately 105 ml of liquid and placed in the ther- 
mostat bath, the latter having previously been 
adjusted in the thermostat bath, the latter 
having previously been adjusted to the required 
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temperature. The exit tubes (J, and J,) of 
the viscometer were connected. 

To expel liquid from bulb (W) to bulb (V) 
(Fig. 4), cocks K, and K, were closed and cock 
K, 
into bulb (W) by opening cock K,. To reverse 
the liquid flow, cocks K, and K, were closed, 
K, was opened, and air pressure was applied to 
bulb (V) by opening K,. The air pressure thus 
applied to (V) was shown on either the mercury 
manometer (M) or the water manometer (N), 
depending on whether cock K, or cock K, was 


was opened. The air pressure was applied 


opened. The water manometer was used, for 
greater accuracy, when the applied air pressure 
was less than 4cm Hg gauge. 

Cock K, was of large bore to enable easy 
exit of the air to the atmosphere from bulb (W) 
as the liquid level in the bulb rose. The various 
cocks and T-joints in the manifold system were 
connected together with short lengths of pressure 
hose which allowed some flexibility and facilitated 
the removal of the viscometer from the thermo- 
With the pressure rubber tubing, 
the viscometer could be adjusted exactly to a 
vertical position in the bath. About 15 minutes 
were required for the sample to attain the 
temperature of the bath, the temperature 
equalizing process being promoted by blowing 
the sample liquid to and fro between bulbs 
(V) and (W). The temperature equilibrium of 
the sample was checked by setting the feed air 
pressure at some arbitrary value and taking 
successive measurements of the time required 
for bulb (V) to be emptied of liquid between the 
graduation levels at (R) and (S). Time measure- 
ments were recorded on a stop-watch reading 
to 0-lsec. Pressure differences could be read 
to + 0-:02cm Hg on the mercury manometer 
and + 0-1 cm water on the water manometer. 

The ordinary Ostwald viscometer has limita- 
tions both as regards the range of data that can 
be obtained from it and also with regard to certain 
assumptions implicit in obtaining the data. 
The modified Ostwald viscometer was developed 
to eliminate these objections and is illustrated 
in Fig. 5. The sample liquid was expelled from 
bulb (V) through the capillary (U) into bulb 
(W), and the time of passage of the liquid level 


stat bath. 
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from the graduation mark (R) to the graduation 
mark (S) was measured. (Z) was a check valve 
with a buoyant float to facilitate the transference 
back of liquid from bulb (W) to bulb (V). 

To make a measurement of viscosity, the 
bulb (V) was completely filled up to cock K, ; 
then with K,, K,, K, and K, closed and K, open, 
the air pressure was applied down the tube (Q). 
In filling the bulb (V), the tube (Q) was also 
partially filled so that with the initial passage 
of air down (Q) displacing this liquid, the liquid 
flow out through (X) into the bulb (W) before the 
liquid in (Q) was displaced beyond the constric- 
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was at one definite pressure, the rate of flow 
calculated from the data was for that pressure 
and was not an average rate of flow derived for 
a pressure range, which was the case when the 
normal Ostwald viscometer was employed. 
Viscometers Nos. I, II, III and IV were used 
to determine the effect of wall slippage, as con- 
sidered in the discussion of the results. The 
solution used was a 0-35°% solution of sodium 
The modified Ostwald 
viscometer was subsequently employed for the 
whole of the viscosity data on the fluids used 
for the experiments on heat transfer [4, 5]. 


carboxymethyl! cellulose. 


tion (R). The time for the liquid to pass through The viscosity data were obtained for the 
700 
Key 
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Fig. 5. Graph of the consistency variables for 0-35% CMC 
solution in water, showing the effect of capillary diameter. 


The data are treated and 
presented graphically in Figs. 6-8 as well as in 
Tables 1 to 17. 


the constriction (R) to constriction (S) 
measured. 

The object in designing the modified Ostwald 
viscometer was to enable the liquid to be discharged 
through the capillary tube under a pressure head 
which was constant and independent of the 
liquid level in either viscometer bulb. It was 
possible with this instrument to take measure- 
ments at pressure heads approaching zero 
without any low limit being imposed by either 
hydrostatic head in the viscometer bulbs. Since 
the whole of the flow for a given measurement 


was following solution. 


(a) 0-25°%, sodium carboxymethyl cellulose 
(CMC) solution ; data at 30°C, 85°C, 42°C, 
49°C, 56°C, 68°C and 70°C. 

(b) 2.5%, polyvinyl alcohol (PVA) solution ; 
data at 30°C, 40°C, 50°C, 60°C and 70°C. 

(c) 0.025%, sodium polymethacrylate (SPM) 
solution; data at 30°C, 40°C, 50°C, 60°C 
and 70°C, 
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Table 1. Viscosity data of 0:25% CMC solution at 30°C. 
Instrument : Modified Ostwald Viscometer. 
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Pressure (cm Hg) Time (sec.) | N’ (cm.) Sg (g/cm sec*) R, (see!) x 10°75 
0-169 2010-2 17-30 6-0 0-147 
0-85 420-6 17°55 33-9 0-69 
1-77 201-5 18-05 68-6 1-45 
242 143-7 18-60 91-1 2-03 
2-84 122-35 19-12 104-0 2-40 
3-44 102-1 19-85 121-2 2-86 
4:39 80-6 20-7 148-4 3-62 
6-04 59-9 22-75 185-6 4-87 
7-09 51-5 23-95 208-4 5-67 
8-29 44-9 25-35 229 | 6-50 

Instrument; Ostwald Viscometer No. 
6-8 19-6 18-8 237-8 | 6-78 
8-05 170 20-1 264-4 7-82 
9-65 14-8 22-5 283-6 | 8-98 
10-6 13-6 23-7 296-1 | 9-77 
11-9 12-35 24-9 316-4 10-77 
12-9 11-5 26-1 827-5 11-56 
14-1 10-7 27-0 345-6 12-41 
14-65 10-35 27-6 351-2 12-84 


Table 2. Viscosity data of 0-25% CMC solution at 35°C. 


Instrument : Modified Ostwald Viscometer. 


Pressure (cm Hg) 


Time (sec) 


1551-0 
253-6 
220-1 
173-7 
129-9 
111-3 
80-6 
70-6 
60-2 
49-8 
44-6 


N’ (em.) 


16-8 
7-55 
18-0 
19-2 
19-35 
20-6 
21-6 
22-75 
25-1 
26-1 


Sg (g/cm sec)* 
10-4 
36-3 
55-4 
70-7 
90°5 
99-8 

134-5 
148-0 
168-0 
193-0 
209-6 


R, (sec) x 10°% 


0-19 
0-83 
1-32 
1-68 
2-25 
2-538 
3-49 
3-99 
4°68 
5-86 
6-45 


235 


| 
/ 
4 
| 
0-25 
0-91 
1-45 
1-87 
2-48 | 
401 
461 
5-46 
6-91 
7-81 | 
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Table 2—{contd.) 
Instrument; Ostwald Viscometer No. IV. 
4-45 28-1 17°55 | 168-0 4-73 
5-3 23-15 | 18-55 | 189-3 5-75 
6-15 20-35 19-2 212 6-53 
6-7 19-0 19-95 222-4 70 
6-8 18-7 20-2 | 223 7-11 
8-0 16-4 21-65 244-5 8-11 
9-15 14-85 23-1 262 8-95 
10-35 13-2 24-7 277 10-06 
11-55 12-2 25-9 295 10-89 
12-1 11-55 26-6 300-5 11-51 
13-3 10-7 28-05 3135 12-42 
13-9 | 10-35 28-55 322 12-85 
| | 
Table 3. Viscosity data of 0-25% CMC solution at 42°C. 
Instrument: Modified Ostwald V iscometer. 

Pressure (cm lig) Time (sec) N’ (em.) Sg (g/cm sec*) R, (sec™') x 1078 
0-48 578-5 17-65 19-0 | 0-50 
0-805 381-8 17-75 31-8 0-76 
1-02 276-1 18-25 39-1 1-055 
1-75 157-5 19-5 62-8 1-85 
2-75 102-7 20-5 94-0 | 2-84 
401 740 22-45 130-1 3-94 
471 64-4 23-5 140-0 | 4-53 
5-41 56-2 25-1 150-7 | 5-19 
5-06 52-05 25-95 160-6 | 5-60 
6-96 45-7 27-7 | 177-0 | 6-37 
7-71 | 41-95 29-9 190-5 6-95 

| 
Instrument ; Ostwald Viscometer No. IV. 
41 26-95 18-85 144-0 4-93 
5-7 20-7 20-35 185-3 6-42 
6-3 18-65 21-35 195-4 7-12 
6-45 18-2 21-55 198-1 7-30 
6-75 17-55 22-1 202 7-57 
7-5 16-0 23-2 214 8-30 
8-25 148 24-45 223 8-98 
9-0 14-05 25-0 238 9-45 
9-6 13-1 26-1 243 10-13 
120 11-05 | 28-75 276 12-02 
12-8 10-5 29-5 285 12-66 
13-4 10-1 | 30-3 292-5 13-16 
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Table 4. Viscosity data of 0-25% CMC solution at 49°C. 
Instrument ; Modified Ostwald Viscometer. 


Pressure (cm Hg) Time (sec) N’ (em.) Sg (g/cm sec*) R, (sec™') x 107% 


0-44 563-4 17-55 17-5 0-52 
0-96 18-6 83-6 1-11 
1-66 19-8 58-7 1-90 
2-46 20-9 79-0 2-75 
2-92 | 21-85 93-5 3-22 
3-32 22-85 102-5 3-67 
4:32 24:3 124-4 4-55 
5-47 | 27-2 140-8 5-59 


Table 5. Viscosity data of 0-25% CMC solution at 56°C. 
Instrument : Modified Ostwald Viscometer. 


Pressure (cm Hg) Time (sec) (em. Sg (g/cm sec*) 
0-34 35- . 13-0 
0-555 399- 21-2 
1-21 43-0 
1-58 . : 54-0 
1-88 . 62-2 
2-28 72-5 
3-18 93-1 
4-23 5- 1142 
528 50: 127-2 
6-23 140-2 


Instrument ; Ostwald Viscometer No. IV. 
42 24-6 20-1 138-3 5-40 
44 | 23-5 20-85 139-4 5-65 . 
5-2 20-8 21-55 159-4 6-39 | 
: 6-05 | 18-35 22-6 177-2 7-24 4 
6-55 17-25 23-15 187-0 7-70 
7-2 | 15-8 | 24-5 194-5 8-41 
7-3 15-6 | 24-65 196-0 8-52 
8-1 | 14-4 25-6 209 9-23 
a+ 9-0 13-2 26-8 222 10-06 
ie 9-45 | 12-85 27-65 225 10-34 
9:85 12-3 28-4 229-5 10-80 
11-0 11-3 29-75 244-5 11-76 
iy 11-65 10-9 30-1 256 12-19 
12-9 9-9 82-2 265-5 13-42 
| 
4 
| R, (sec) x 107 
0-46 
0-73 
1-55 
2-00 
2-38 
2-88 
3-85 
4-78 
6-45 
| | | 
287 
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Table 5—{contd.) 
Instrument ; Ostwald Viscometer No. IV. a 
| 7 
3-6 26-0 21-05 118-0 5-11 7 
45 21:8 22-6 181-8 6-18 
5-25 | 19-0 28-95 145-0 00 a 
6-8 16-6 25-65 162-5 8-00 
ay 7-8 14-95 26-6 181-4 8-89 
7-6 14:3 27-7 1810 9-30 
8-5 13-1 29-25 191-8 10-14 
9-65 12-0 30-45 209-6 11-08 
10-55 11-05 82-4 215-7 12-08 
12-2 9-8 34-3 235-5 13-56 


Table 6. Viscosity data of 025%, CMC solution at 63°C. 
Instrument - Modified Ostwald Viscometer. 


Pressure (cm Hg) Time (sec) N’ (em.) Sg (g/cm sec*) R, (sec™') x 10°75 
0-24 652-2 18°3 9-2 0-38 
415-6 18-3 16-8 0-70 
0-75 19-55 26-8 1-12 
1-41 138-5 21-75 45-4 2-11 
2:18 95-1 23.5 65-0 8-07 
3-04 70-9 26-8 79-5 411 
3-34 66-8 28-2 83-0 4°37 


4°69 31-3 104-9 5-66 


Instrument Ostwald Viscometer No. 


3-55 24-2 22-8 102-8 5-49 


3-95 22-45 23-45 111-2 5-91 
42 21-3 24-0 115-8 6-24 
4-55 20-3 24-6 122-2 6-65 
5-2 18-15 25-85 182-9 7-32 
5-9 16-3 28-0 139-1 8-15 
6-2 15-15 28-2 145-2 8-38 
70 14:3 30-15 153-6 | 9-29 
7-45 13-9 30-1 163-7 9-56 
8-15 12-95 31-2 172-6 10-26 
85 12-4 32-6 172-7 | 10-72 
9-1 11-8 33-2 181-2 11-26 
10-1 10-85 34-7 192-5 12-25 
11-05 10-18 35-6 205 | 18-11 
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Table 7. Viscosity data of 0-25% CMC solution at 70°C. 
Instrument : Modified Ostwald V iscometer. 


: Pressure (cm Hg) Time (sec) N’ (cm) | Sg (g/em sec®) | R, (sec) x 1078 

0-32 489-5 20-5 10-9 | 0-60 
0-70 231-7 211 23-2 1-26 

1-41 124-7 23-1 42-7 2-34 

1-67 108-7 23-8 49-1 2-68 

2-50 77-9 27-2 60-1 8-75 

3-10 65-3 29-9 72-6 4-47 

3-75 56-1 32-8 80-0 5-20 

4-30 50-8 33-3 90-5 5-74 

Instrument ; Ostwald Viscometer No. IV. 
| 

2-35 32-0 21-7 71-6 | 4-15 

415 20-15 26-9 102-7 | 6-60 

485 18-25 27-65 115-8 7-28 

5-4 16-85 28-8 124-0 7-89 

59 155 30-35 128-4 | 8-57 

6-85 13-95 31-95 141-9 9-53 

7-2 13-3 33-4 | 142-6 9-99 

7-8 12-55 34-6 149-1 10-59 

9-1 11-18 36-5 165-0 11-92 

10-5 9-95 39-0 178-0 13-36 

Table 8. Viscosity data of 2-5% PV A solution at 30°C. 
Instrument ; Modified Oswald Viscometer. 
Pressure (cm Hg) | Time (sec) N’ (em.) Sg (g/cm sec™*) | R, (sec™') x 107% 

0-71 415-0 | 17-6 28-2 0-70 

1-17 254-6 17-95 45-6 1-14 

1-46 206-2 18-30 55-9 1-41 

2-02 | 156-1 18-95 74-6 1-87 

2-84 114-0 19-6 101-4 2-56 

3-34 96-6 20-2 115-7 3-02 

3-89 84-6 20-75 131-2 8-45 

4-49 75-4 21-35 147-3 3-87 

5-15 66-3 22-15 162-8 4-40 

6-04 58-6 22-9 184-7 4-98 

6-14 57-8 23-0 187-0 5-05 

7:14 50-85 24-25 206 5-74 

8-04 46-7 24-8 227 6-25 

8-59 44-0 25-85 287-5 6-68 


q 


Ju Cum Cuvu, K. G. Buraiver and Frank Brown : The viscosity of pseudo-plastic fluids 


Table 9. Viscosity data of 2.5% PV A solution at 40°C. 
Instrument : Modified Ostwald V iscometer. 


Pressure (cm Hg) Time (sec) | R, (see) x 10% 


0-59 398-1 0-73 
0-74 316-3 . . 0-92 
1-16 207-9 1-41 
1-64 149-7 1-05 
2-35 108-05 . 2-70 

89-2 . 3-27 
3-4 78:3 3-72 
3-8 71-9 23-25 4-06 

63:1 241 126-3 4-62 
47 60-7 24-7 133-1 4-80 
5-65 52-4 26-6 148-6 5-57 

48°35 27-8 157-4 6-04 
6-6 47-2 27°5 168-0 6-17 
72 43-9 29-0 1740 6-64 


Table 10. Viscosity data of 25%, PV A solution at 50 
Instrument : Modified Ostwald V iscometer. 


Pressure (cm Hg) Time (sec) | j Sg (g/cm 
434°8 16-3 
268-0 25-9 
175-4 37-4 
126-5 51-4 
103-0 60-9 

92-6 64-6 
89-5 68-3 
78-55 75-8 
69-4 86-1 
104-0 
51-9 113-2 
47°5 117-9 
45:1 125-1 
42-7 34-9 126-4 


Table 11. Viscosity data of 2.5% PV A solution at 60 
Instrument Modified Ostwald Viscometer. 


Time (sec) N’ (em) Sg (g/em 


Pressure (cm Hg) 


622-1 19-85 8-8 
261-3 20-5 
22-25 
122-0 23°5 
93-5 25-05 
83-6 25-8 
72-05 25-4 
66-6 20-5 
59-1 32-4 
52-95 83-5 
45°8 36-9 
42-15 387-6 
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a 
= 3 
1954 
sec?) R, x 10° 
0-67 
1-09 
| 
3-05 
3-26 
3-71 
one 
6-14 
4 6-47 
6-83 
| 
R, (sec") x 10% 
0-25 0-47 
0-64 1-12 
1-05 1-86 
| 2-39 a 
201 8-12 q 
2-36 3-49 
2-71 4-00 
3-08 | 4-28 
3-61 4-04 
+21 5-51 
5-26 6-87 
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Table 12. Viscosity data of 2.5% PV A solution at 70°C. 

Pressure (cm Hg) Time (sec) N’ (cm) Sg (g/cm sec*) R, (sec) x 10% 
0-28 492-4 21-2 9-2 0-59 
0-43 316-3 218 13-8 0-92 
0-94 156-2 24-0 27-5 1-87 
1-23 123-9 25-3 34:1 2-36 
1-41 109-5 27-2 36-2 2-66 
1-91 89-05 28-4 47-1 3-28 
2-13 81-6 30-0 49-7 3-58 
2-59 72-0 31-1 58-8 | 4-05 
3-26 60-8 34-0 66-8 4-80 
3-46 58-3 34-3 70-6 501 
$31 50-5 37-4 80-5 5-78 
1-6 38-1 85-8 6-08 
501 415-3 40-7 86-1 6-42 


Table 13. 
Instrument : Modified Ostwald Viscometer. 


Viscosity data of 0-02 


SPM solution at 30°C. 


Pressure (cm Time (sec) N’ (em) Sg (g/cm sec?) (sec) 
O-87 422-4 16-65 36-6 0-69 
1-19 290-9 17-2 48-4 1-00 
1-47 230-5 17-6 58-4 1-27 
2-06 153-8 18-75 76-9 1-90 
2-42 130-6 19-3 87-8 2-24 
2-92 101-9 20-4 100-2 2-86 
3-52 83-05 21°5 114-6 3°51 
4-07 708 22-9 124°8 4°12 
4-42 65-6 23-6 131-1 | 4°45 
5-02 58-7 24-7 142-2 4-97 
5-27 56-5 25-3 145-7 5°25 
6-22 48-1 27-7 157-1 6-07 
6-72 45-0 29-3 160-5 6-48 
7-37 41-9 30-4 160-1 6-96 

Table 14. Wiscosity data of 0-025% SPM solution at 40°C. 


Instrument : Modified Ostwald Viscometer. 


Pressure (cm Hg) Time (sec) N’ (em) Sg (g/cm sec*) R, (sec) x 10% 
0-64 465-1 17-3 25-9 0-63 
0-80 365-9 17-7 31-6 0-80 
1-36 211-1 18-55 51-3 1-38 
1-65 161-1 19-3 60-1 | 1-81 
2-13 123-6 20-3 73-1 2-36 
2-63 98-2 21-5 85-5 2-97 
3-08 82-6 22-7 94-9 8-53 
3-48 73-7 23-5 103-6 3-96 
3-98 64-7 24-7 112-8 4°51 
4-73 56-1 27-2 121-8 5-20 
5-87 50-7 28-8 130-7 5-75 
5-82 47-7 29-9 136-4 | 6-11 
6-43 43-6 32-8 143-0 6-69 
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Table 15. Viscosity data of 0-025% SPM solution at 50°C. 
Instrument : Modified Ostwald Viscometer. 
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N’ (em) 


18-25 
19-3 
20-3 
20-9 
22-6 
24-2 
25-4 
26-4 
27-1 
29-3 
30:8 
33-4 
33-5 


| Sg (g/cm sec*) 


22-3 
38-1 
49-0 
57-6 
70-7 
80-4 
87-6 
93-5 
99-5 
103-7 
112-1 
117-0 
118-6 


Table 16. 


Viscosity data 


Instrument : Modified Ostwald Viscometer. 


of 0-02 


~/o 


SPM solution at 60°C. 


Pressure (cm Hg) 
0-51 
0-75 
0-98 
1-33 
1-85 
2-23 
2-79 
2-89 
3-19 
3-66 
4-29 
504 


Time (sec) 


359-7 
248-7 
190-0 
145°1 
116-0 
102-2 

86-4 

718 

710 

65°5 
57-6 
45-0 


N 


(cm) 


19-5 
20-0 
20-7 
21-65 
23-1 
24-0 
25-7 
28-3 
28-0 
29-6 
31-6 
85-0 
37-3 
37-9 


Sg (g/cm sec*) 


18-3 
26-2 
33-1 
43-0 
49-6 
53-8 
60-7 
69-0 
72-3 
75-5 
80-9 
85-9 
94-6 
97-6 


R, (sec") x 10% 


0-83 
1-17 
1-54 
201 
2-52 
2-85 
3-38 
4-06 
411 
4°46 
4-08 
5-66 
6-27 


6-48 


Table 17. 


Viscosity data of 0-025°%, SPM solution at 70°C. 


Instrument : Modified Ostwald Viscometer. 


Pressure (cm Hg) 


0-24 
0-47 
0-74 
102 
1-26 
1-72 
2-24 
234 
2-54 
2-80 
3-29 
3-49 
3-69 
8-99 
474 


Time (sec) 


505-3 
303-9 
197-5 
1445 
117-7 
99-3 
76-85 
756 
69-6 
65-2 
56-9 
55°5 
49-8 
448 


N 


(cm) 


20-7 
21-7 
22-8 
24°4 
26-9 
27:8 
32-4 
32-3 
33-8 
343 
38-0 
38-3 
42-3 
43-9 


Sg (g/cm sec?) 


8-1 
15-2 
22-7 
29-2 
32:8 
43-3 


R, (sec 10° 


0-49 
0-906 
1-48 
2-02 
2-48 
3-04 
3-80 
3-86 
419 
4-48 
5-07 
5-26 
5°45 
5-85 
6°51 
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q 
(com Hg) (sec) | R, (sec) x 10% 
1-05 219-9 | 
1-42 160-4 
172 133-5 | 182 
2-28 97-6 
810 3-60 q 
3-18 72-4 00 
8-58 66-05 
4-33 55-8 | 
4-03 50-8 5-74 q 
5-58 46-1 6-34 
: 
| 1954 
| 
5-29 | 
| 
= : 
| 
| 48-4 
| 
52-6 
| 
60-6 
63-6 
66-6 
| 66-0 | | 
| || 


In applying the Hacen-PoisevILLe equation 
to the results obtainable from an Ostwald type 
viscometer, certain sources of error arise, these 
were discussed in some detail by Haut and Fuoss 
[7]. They considered the deviations to arise 
as follows : 


1. All the work done by the driving pressure 
is not dissipated as heat in overcoming the 
viscous friction. 


2. Flow near the ends of the capillary tube 
is not a laminar. 

8. The volume of liquid which flows through 
the tube is not the volume of the dry bulb. 


A method for the analysis of data for a pseudo- 
plastic fluid to eliminate these errors was 
considered for the present work, but the results 
do not lend themselves to extrapolation in the 
manner recommended. Hence the following 
method of calibration was employed : 


1. The volume of the viscometer bulb (V) 
between the constrictions at (R) and (S) 
was measured by filling the bulb with water 


360) 


+ 


| | 


T 
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Fig. 6. Temperature effect upon shear rate and stress relation 
for 0-25 CMC solution in water, with temperature as parameter. 
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draining off the water, and weighing it. 
The mean of several determinations was 
used in subsequent calculations. The nipple 
(Y) on Fig. 4 facilitated this operation. 

2. The diameter of the capillary tube was 
measured, prior to the construction of the 
viscometer, by two methods : firstly, by 
direct measurement of the diameter at 
either end with a travelling microscope ; 
secondly, by measuring the length of a 
weighed thread of mercury at various points 
along the bore of the capillary. 


8. The effective or equivalent length (N’ cm) 
was obtained by performing viscosity 
measurements with solutions of pure cane 
sugar, the viscosity of which was known (9). 
The equivalent length was then calculated 
from the relationship. 


(9) 


In practice, the equivalent length was found 
to vary with both rate of flow and with viscosity 


rT 


| 
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to vary with both rate of flow and with the 
viscosity of the fluid. Thus it was necessary to 
established the equivalent length over a range of 
values for both of these variables. The full 
procedure is described elsewhere (4). 

The data actually recorded were the pressure 
applied to bulb (V) as shown on the manometer 
(M) or (N) in Fig. 3, and the time required to 
empty the bulb (V). From these a plot of applied 
pressure versus R, was made (4). The applied 
pressure was expressed in cm Hg. In 
of the modified Ostwald viscometer an amount 


the case 


graph of N’ versus R, the required value of N’ 
was that which gave the same value for r, 
for the datum value of R,, as that just estimated on 
the graph of P versus RK, With N’ known, 
the shearing stress at the wall, Sg, for any given 
viscosity determination could then be calculated 
from the relation 


Poy 


Sg (10) 


Discussion of results 


As the precautions required to ensure that the 


equivalent to the hydrostatic pressure head from viscometers were accurately calibrated were 
240 T T 
Key 
DaTA aT 3O°C © 
200 40°C x = 
© 
60°C x 
70°C © 
EEO 4 
Sq 
| | 
190 
| 
80 
40 


Fig. 
for 2-5°% PVA solution in water, with temperatureas parameter. 


(R) to the level at (X) in Fig. 4 was subtracted 
from the manometer reading to give the effective 
pressure drop across the viscometer. To calculate 
the p 
fluids the procedure was as follows: R, 
the 
conjunction with the known values of V and R. 


resulting viscosities with pseudo-plastic 


was 
calculated from measured values of p» in 
A graph of P versus R, was constructed from data 
given [4, 6] and on this graph the point corres- 
ponding to the data values of P and R, was 
located and the value of the parameter + corres- 


ponding to this location was estimated. From 


Temperature effect upon shear rate and stress relation 


somewhat elaborate, it is appropriate to mention 
briefly the theoretical background which necessi 
tated this. 

It has been fairly well established that the 


7 slip * occurs in a flow of pseudo plastic fluid 
in a capillary tube. The term “slip” was des- 
cribed by Winpinc et al. [7, 14, 16]: “For a 


the 
until a 

the 
independent of diameter. This type of behaviour 


(16]"". The 


given rate of shear, shearing stress is a 


function of diameter critical value is 


reached, above which shearing stress is 


has been designated as ‘slip’ 
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viscometers of different diameters are used to 
determine the relationship between the shear 
rate and shearing stress. It was found that the 
minimum diameter of the capillary tube was 
0-199 cm in order to be free from the “ slip.” 
The viscometer of this minimum diameter was 
used in the determination of viscosity of the 
pseudo-plastic fluids. 

Baker [2] established that turbulence losses 
due to enlargement, contractions, and orifices 
for pseudo-plastic fluids substantially 
similar to the corresponding losses for the simple 


were 


errors due to slip were present. Thus a certain 
margin of safety was introduced into the choice 
of the minimum suitable diameter. 

Fig. 6 shows that providing R, was less than 
7000 sec’ the capillary of 0-199cm diameter 
met the required condition. In the paper on heat 
transfer (5) it will be seen that in the application of 
the viscosity data to the heat transfer correlations 
the requirement that R, should be less than 
7000 sec™' was always observed. 

Referring to Figs. 6, 7 and 8, it is seen that the 
range of data are for R, = 0-7000 sec for the 
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Fig. 8. 


0-025% SPM solution in water, with temperature as 


Newtonian fluids. 
effect of the capillary tube, the viscometer of 
minimum diameter was calibrated with a simple 
Newtonian fluid. 

Fig. 6 shows the result of a test to establish 
the minimum required diameter. 0.35°,, CMC 
This solution 


To compensate for the end 


solution was used for this test. 
was of a somewhat higher concentration of CMC 
than was used in later experiments. It was so 
selected to give a viscocity higher than that 
anticipated in the remainder of the work, in 
order to provide a more severe test of whether 


Temperature effect upon shear rate and stress relation for 


parameter. 


PVA and SPM solutions, and R, = 0-14,000 sec 
for the CMC solution. The data in all three cases 
for R, = 0-7,000 sec"! were obtained on the modi- 
fied Ostwald viscometer, but the additional data on 
the CMC were obtained on Viscometer No. IV. 
The object of taking the additional data on 
CMC was two-fold. Firstly, it provided a test of 
whether the shearing rate vs. stress curve derived 
from data on Viscometer No. IV was a continua- 
tion of that derived from the modified Ostwald 
viscometer, thus giving a check on the satisfactory 
operation of this design of instrument. Secondly, 
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it was desired to have a fairly wide range of data 
in an attempt to obtain graphically the limit, 
#..» to which the slope of the shearing rate vs. 
stress curve approaches as Sg becomes infinity. 
The results of this attempt were not too satis- 
factory, and consequently the use of ». in 
correlating the heat transfer data [5] was avoided. 
In view of this conclusion the viscosity data for 
the PVA and SPM solutions were only taken 
over a more limited range. 

It was preferred to established an analytical 
relation between the shear rate and the stress 
at different temperature, rather than to read 
the date necessary for the heat transfer calcula- 
tions directly from graphs. By so doing, more 
accurate values could be obtained for low values 
of R,. 
the consistency 
intermediate temperatures. 

None of the standard methods of rectification 


It also facilitated interpolation between 


curves for data required at 


appeared to give satisfactory results over the 
whole range of the viscosity data, although a 
log-log plot such as that used by Krircer and 
Manon [9] resulted in a straight line plot over a 
limited range. It was thus necessary to consider 
an empirical equation. The hyperbolic form 
used previously by WILLIAMSON [15] and WINDING 
et al. [16] was found satisfactory over a sufficient 
range for the present work. 

Considering first the CMC data, it is seen that 
the dotted lines through the data points at higher 
Hence 


" may be determined graphi 


values of R, are substantially straight 
SoZ 
cally. Use of the values thus obtained, plus 


lve for the 


and 


one other datum point, suffices to so 
constant “a” and to give an equation of the 
form used by WINDING. 

however, did not 


The resulting equation, 


satisfactorily predict the data values for the 
low values of R.. An alternative method of 
fitting the WinDING equation to three selected 
point values spaced over the range 0 R. 
7,000 sec.' gave a very good prediction of the 
data over the whole of this range, although the 
resulting curve inaccurate values for 
R, >» 7.000 sec 


This deviation may be due to three 


gave 


possible 


causes 
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1. Viscosity determinations at high values 
of R, were difficult to measure accurately 
due to the very short time required to 
empty the viscometer bulb. 

2. Due to the “slip” effect, the effect of 
capillary diameter on the slope of the 
shearing rate vs. stress curve is no longer 
negligible at high values of R,. 

3. WunpiNnc’s equation is empirical and there 
may be no theoretical reason why it 

should hold over the whole range of R,. 

The so called WinpiNG equation is of the same 
form as that originally proposed by WILLIAMSON as 
given by eq. (7). Since the actual shearing stress 
and shear rate between parallel plates was assumed 
in Wituiamson’s original method, the rate of 
shear R. calculated at the wall as suggested by 
Winpinc is different from the actual one between 
the parallel plates. Any method using apparent 
viscosity alone is bound to be less accurate than 
those using differential viscosity, though all 
previous authors have used the former as the 
easy way out. 

The analytical expressions to fit the CMC 
solution viscosity data at the various tempera- 
tures having been derived, it was graphically 
possible to read off the values of the constants 
for the Winding equation at any given tempera- 
ture and so compute the apparent viscosity 
for the appropriate value of R. The same 
procedure was followed in the development of the 
viscosity data for the PVA and SPM solutions. 

It should be 


approach was selected solely because it represented 


emphasized that the above 


a more convenient and accurate method of 


obtaining the viscosity data required for the 
heat transfer correlation. The WinpDING equation 
is empirical, but its application in the present 
work is justified by the fact that it was only 
used over the range of data for which it was 


shown to hold. 


ACKNOWLEDGEMENTS 


Grateful Prof. 
Winpinc of Cornell University, and Prof. Orumer 


acknowledgement is made _ to 
«f Polytechnic Institute of Brooklyn, for their 
comments and interest ; to ECA for the Exchange 
Fellowship awarded to two of the authors. 


46 


VOL. 
3 
1954 
- 


PHeg 


NOMENCLATURE 


Definition 
Diameter of tube 
Mass velocity 
Gravitation constant 
Length of tube 


Equivalent length of capillary 
tube 


Pressure differential along tube 
Inside radius of tube 
Reynolds number 


Apparent Reynolds number 
Rate of sl 
e of shear | — 


Shear stress at wall 


A constant in the WILLIAMSON 
and WINDING equations 


Volume rate of flow 


Linear velocity along the tube 
Density of mercury 
Time of efflux 


Coefficient of Viscosity of New- 
tonian fluid 


Coellicient of apparent viscosity 


Coeflicient of differential vis- 


cosity 


constant in the 


and WINDING equations 


Units 

em. 

g/sec. (em.)? 
g/sec. (em.)* 


em. 


em, 


em of mercury 


em 
dimensionless 
sec.~! 


g—mass /(sec)” 


g/em. sec.” 


As indicated 
in text. 


ml /second 
g/ml 


secs, 


/ (em) (sec) 


same 


same 


same 
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The mechanism of mass transfer of solutes across liquid-liquid 
interfaces 


Part I: The determination of individual transfer coefficients for binary systems 


J. B. Lewis 
Chemical Engineering Division, Atomic Energy Research Establishment, Harwell, Didcot, Berks. 


( Received September 1954) 


Summary—aA new type of transfer cell has been designed and used for the determination of 
the individual mass transfer coefficients, by the Colburn-Welsh technique, for eight partially 
miscible binary systems, in six of which one of the phases was water. It was found that the 
transfer of mass to the interface took place by eddy not molecular diffusion. The individual 
coefficients for all phases, except furfural, were correlated on one curve by plotting , /v,; against 
Re, + Reg n2/4,. This curve could be represented with a slight loss of accuracy by the expression : 


The effect of adsorbed interfacial films was examined and it was found that while mobile films 
had little effect rigid protein films caused a retardation of transfer, probably by damping inter- 
facial turbulence. 

Since the transfer coefficients are dependent only on the Reynolds numbers and physical 
properties of the two phases, the correlation enables individual coefficients to be calculated 
for immiscible solvent-water pairs where the mutual saturation technique cannot be employed. 


Résumé—Un nouveau modeéle de cellule de transfert a été congu et utilisé dans la détermination 
du coefficient individuel du transfert de masses par le procédé CoLpurN-Wetsn dans le cas de 
huit systémes binaires partiellement miscibles pour six desquels l'eau était une des phases. On 
a trouvé que ce transfert de masse a l'interface se propduit par tourbillon et non par diffusion 
moléculaire. 

Pour toutes les phases, excepté le furfural, les coefficients spécifiques ont été portés sur la 
courbe k, /v,; + (Re, + Reg Cette courbe pourrait étre représentée avec un peu moins 
de précision par l'expression : 


60k, /v, + O76. (Re, + Reg 


On examine leffect de couches absorbées & l'interface et quoique les couches mobiles aient 
un faible effet, on trouve que les couches rigides de protéines provoquent un ralentissement du 
transfert probablement en atténuant la turbulence a l'interface. 

Puisque les coefficients de transfert dépendent seulement des nombres de Reynolds et des 
propriétés physiques des deux phases, la formule permet de calculer les coefficients pour des 
binaires eau-solvant immiscible dans le cas ot: le procédé de saturation réciproque ne peut étre 
utilisé. 


The transfer. of solutes, particularly organic transfer is followed by a homogeneous chemical 
solutes across liquid-liquid interfaces is usually reaction the reaction is considered to take place 
considered to be diffusion controlled. Normally at a zone inside the relevant film. This two film 
both phases are in turbulent motion but it is theory, adapted by Wurrman [1] from the general 
assumed that films or layers exist in each phase boundary layer theory of Nernst [2], is used 
immediately adjacent to the interface where the extensively to correlate the performance of 
flow is laminar, and that diffusion through these extraction equipment. Although various aspects 
layers is slow and rate-determining. When the of the theory have been criticised no experiments 
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have been carried out to test its validity directly. 

Van Name [3] and others have shown that when 
a chemical reaction occurs at a solid-liquid inter- 
face then the process will be chemically controlled 
if the reaction is very slow, diffusion controlled if 
the reaction is very fast and of an intermediate 
complex nature if the chemical reaction rate is 
of the same order as the diffusion rate. A similar 
state of affairs must hold for reactions at liquid- 
liquid interfaces. Since diffusion and first-order 
chemical reactions have similar kinetics a slow 
first-order chemical step will be equivalent to 
diffusion through a fictitious film. In such systems 
therefore, the two-film theory may be used to 
correlate data provided that the limitations of 
the theory are remembered. In general, however, 
the two film theory should only be used where 
transport to and from the interface is_ the 
controlling factor. 

In recent years the hydrodynamics of laminar 
films at solid-liquid interfaces have been worked 
out rigorously [4]. The behaviour of films at 
gas-liquid and liquid-liquid interfaces is less well 
known but there are indications that turbulence 
may in some cases persist up to the interface 
itself [5]. Transport of matter to the interface 
will then be at least partly by eddy diffusion. 
Since eddy diffusion is much faster than molecular 
diffusion there is thus more chance of any chemical 
step being rate-determining. Chemical steps are 
known to occur in the extraction of many inor- 
ganic solutes and DrickamMer has shown that 
they may even be involved in some apparently 
simple organic systems [6]. Information as to 
the mechanism of transfer to the interface and 
on the magnitude of any chemical step is clearly 
needed if the performances of different types of 
extraction equipment are to be compared. 

No work has been reported on the nature of 
transport to liquid-liquid interfaces but a few 
papers have been published on the extraction of 
inorganic and organic solutes from water into 
organic solvents [6, 7, 8]. In these experiments 
stirring was either absent or else an attempt was 
made to standardise the effect by making obser- 
vations at a number of stirring speeds and extra- 
polating back to zero stirring rate. When stirring 
was absent effects due to local heating, density 


gradients, etc., became very important, and the 
results were sometimes suspect. In some experi- 
ments local stirring was prevented by packing 
the bulk phases with wire gauze [6]. Separate 
experiments were then necessary to determine 
the unknown diffusion path lengths. The stirred 
experiments [7] in general have suffered from non- 
reproducibility of the turbulence, insufficient 
attention having been paid to the hydrodynamics 
of the apparatus. The method adopted of extra- 
polating the stirring rate back to a fictitious zero 
rate has recently been shown to be at fault [8]. 

In the present work the unsteady-state transfer 
of solutes between phases was investigated using 
a specially designed transfer cell in which the 
degree of turbulence could be accurately con- 
trolled. Initially the rate of mutual saturation of 
a number of partially miscible binary systems was 
measured so that the individual transfer coeffi- 
cients for these phases could be obtained by the 
CoLBuRN and method [9]. The depend- 
ence of the individual coefficients on the Reynolds 
numbers of the two phases, on the diffusion 
coefficients and other physical properties of the 
two components and on the temperature, was 
examined and a method is derived for the approxi- 
mate correlation of the data. The results of this 
investigation are summarised below. In Part 
II [10], this work has been extended to consider 
the rates of transfer of organic solutes between 
solvent and water phases. 


EXPERIMENTAL 

Apparatus 

The experiments were carried out with a cylindrical 
transfer cell in which each phase was independently 
stirred by a central paddle. Preliminary experi- 
ments were first carried out in simple glass vessels 
to determine the optimum conditions for pro- 
ducing reasonably smooth interfaces of repro- 
ducible area at high stirrer speeds. The vessels 
were closed at the top and completely filled in 
order to prevent cavitation of the upper liquid 
surface. It was found necessary to restrict the 
interface to an annular gap between a central 
circular baffle and a circumferential wall baffle. 
the central baffle preventing cavitation and the 
outer baffle reducing irregular wall effects to a 
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minimum. The two halves of the cell were made 


geometrically similar so that the film coefficients 
for the two phases could be correlated together. 
Various sizes and shapes of stirrer paddles were 
tried and the dimensions of the baffles were 
varied until a satisfactory stirring range could be 
obtained without break-up of the interface. Both 
baffles were bevelled to enable drops of either 
phase accidentally introduced into the other to 
roll to the interface. 


c 
L 
A FILLING PLUG 
B POLYTHENE PLUG SUPPORTING ELECTRODE 
STIMRER 
D 
E wTeRFace 


Fig. 1. Transfer cell. A — filling plug, B — polythene plug 
supporting electrodes, C — stirrers, D — baffles, E 
face. 


inter- 


The final cell, shown in Fig. 1, consisted 
essentially of a cylinder 4in. in diameter and 
8 in. high constructed of two 1} in. lengths of 
4 in. i.d. glass pipe clamped between two flat end 
plates. These and all other metal parts inside 
the cell were made of stainless steel. The two 
glass sections were separated by the circum- 
ferential baffle, and this together with the central 
baffle divided the cell into two identical halves 
each of ca 300 ml. volume, the annular gap area 
being about 30sq.cm. Gaskets suitable for the 
system being examined, e.g. rubber, P.V.C., and 
polythene-sheathed rubber, were inserted between 
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the various parts. The bottom plate was solid 
but the upper had three circular holes, a central 
hole for the stirrer shafts, a filling hole normally 
closed by a plug and a hole accommodating a 
polythene plug supporting two platinum elec- 
trodes. Both plugs were made to finish flush with 
the inner surface of the upper plate, and two 
small holes in the filling plug enabled hypodermic 
needles to be inserted to sample the two phases. 
The electrodes were connected to a specially 
constructed conductivity bridge enabling the 
conductivity, and hence the concentrations, of 
suitable inorganic solutions to be determined 
during an experiment. To reduce interference 
with the turbulence in the upper compartment 
the electrodes were parallel to the walk and 
curved concentrically with the sides. 

Three concentric shafts passed through the 
centre hole, an inner rotating shaft on to which 
was fixed the stirrer for the lower compartment 
a fixed central shaft supporting the central baffle 
at the interface, and an outer rotating shaft 
driving the upper compartment stirrer. The two 
stirrers were driven by } h.p. D.C. motors con- 
trolled by commercial electronic speed controllers 
which enabled the speeds of both stirrers to be 
varied independently from 0 to 1,400 r.p.m. 
Reversing switches were included in the circuit. 
The cell was suspended in a water thermostat to 
enable experiments to be carried out at various 
accurately-controlled temperatures. 


Systems used 


Mutual saturation runs were carried out with 
six different solvent-water systems and also with 
two non-aqueous systems. For the solvent-water 
experiments the solvents were aniline, ethyl 
ethyl aceto-acetate, ethyl formate, 
furfural and iso-butanol, while formic acid-benzene 
and aniline-cyclohexane were used as the non- 
aqueous systems. The best available grades of 
solvents were used and with the exception of 
formic acid, which was purified by fractional 
freezing, the solvents were all redistilled in a 
1 in. i.d. laboratory column packed with 8 ft of 
6 mm Raschig rings. Furfural was distilled under 
filter pump vacuum. 

The physical properties of the compounds used 


acetate, 


5 i” 
Nels 
tle 
1954 
rey 


J. B. Lewis : The mechanism of mass transfer of solutes across liquid-liquid interfaces — I 


are tabulated in Table 1, much of the data having 
been determined experimentally. No diffusion 
coefficient data were available for the solvent- 
water systems, and most of these were therefore 
determined experimentally (full details of this 
work will be published elsewhere). The dimension- 
less Schmidt groups (» pJ)) were calculated from 
Table 1. Inter- 
facial tensions of the actual phases used were 


these values and are included in 


determined by the pendant drop method [11]. 


Procedure 
Before assembling the cell for an experiment each 
part was carefully cleaned with a suitable solvent, 


rinsed in redistilled acetone and dried in a hot air 


steeped in chromic acid. It was necessary with 
some systems to soak the P.V.C. or rubber gaskets 
in the clean solvent beforehand. Great care was 
taken to avoid contaminating the cell with grease 
or surface active agents. The cell was assembled 
immediately before use, replaced in the bath and 


had 


brought to the correct temperature in a thermo- 


the two phases, which previously been 
stat bath, added in such a manner that neither 
The 


heavier phase was run in until the top surface 


phase wetted the wrong compartment. 


was somewhat above the mid-point of the annulus 
so that on taking the initial sample from this 
phase the interface fell to the correct level. 
Sufficient of the lighter phase was then run into 


oven. Every few days the components were the top compartment to fill it completely and to 
Table 1. Physical properties of solvents. 
Schmidt Group 
Te mp. p n 
Solvent gm In water Water in dynes /cm Remarks 

Furfural 20 1-16 1-70 960 1,600 6 

1-15 1-40 620 720 

Ww 1-14 1-23 - } 5 

50 1-13 1-11 
Ethyl acetate 20 0-90 0-46 1,000 160 9 

30 0-89 0-40 +60 130 - 

0-88 0-37 - 6-5 

50 0-86 0-35 
Ethyl formate 20 0-92 0-40 1,000° 100* 10-5 *Estimated 

values 

Ethyl aceto-acetate 20 1-03 1-71 | 1,400° 1,800* 4-7 
Aniline 20 1-02 1,090 6,200 5-8 
Isobutanol 20 0-80 3-95 | 1,200 14,000 2-1 

50 0-78 1-65 . - 

| 
Formic acid 20 1-21 1-80 - | - ~ 
Benzene 20 0-88 0-65 - 
| 

Cyclohexane 20 0-80 0-96 - - - 
Water 20 1-00 1-00 

80 1-00 | 0-80 

Ww 0-99 0-66 

50 0-99 0-55 
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rise part-way up the filling and centre holes. On 
inserting the filling plug the excess volume of 
about 2 c.c. was forced up into the centre hole. 
This excess volume acted as a reservoir when 
samples were removed from the upper compart- 
ment so that the level did not fall below the 
top plate. Although this excess liquid was not 
adequately stirred the error introduced by its 
somewhat different composition was negligible. 
Due to the variable thickness of the gaskets the 
volume of the two phases varied by a maximum 
of + 2c.c. 

The stirring rates necessary to give the correct 
Reynolds numbers had been previously calculated 
and the controllers set at this speed. During each 
run the motor speeds were checked by a tacho- 
meter and frequently the stirrer rates were 
further checked by a stroboscope. 

Both phases were sampled immediately after 
starting the stirrers and subsequently at intervals 
of 3-5 minutes, by inserting long hypodermic 
needles through the filling plug and drawing 
samples from the middle of each phase. Large 
samples of each phase were also taken at the end 
of each run. All samples but the last were kept 
as small as possible consistent with accuracy of 
analysis, but when only drop samples were 
required for refractometric analysis at least 1 ml. 
was taken into a syringe, and after shaking, three 
or more drops were taken from this. The remain- 
der of the sample was then returned into the 
cell. If larger samples were required, as e.g. for 
Kari Fischer titrations of water content, these 
were kept to a minimum during the run, only 
three or four intermediate samples being taken, 
and synthetic solutions of approximately the same 
composition were injected so as to keep the top 
compartment full and the interface level con- 
stant. The runs normally lasted about 40 minutes. 

A number of transfer experiments were carried 
out with the viscosity of the solvent phase 
increased by the addition of collodion. The 
viscosities were determined before and after each 
run. Transfer runs were also carried out with 
both mobile and rigid adsorbed films of surface- 
active agents present at the interface. The cell 
was filled in the usual manner, the run started 
then the stirrer stopped and the surface active 
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agent, dissolved in a suitable solvent, injected on 
to the interface by means of an “Agla” precision 
syringe connected to a long needle. A long chain 
quaternary ammonium compound C,,H,,N(CH,),- 
Cl was used for the mobile film, this being 
insoluble in either phase.* Enough of a petroleum 
ether-propanol solution was added to give either 
1 or 4 molecules per 150A? of surface area. For 
rigid films excess of a protein, gliadin, in ethyl 
alcohol solution was used. In all cases the 
volumes of the solutions added were negligible. 


Analysis 


Analysis of both phases in the aniline-cyclohexane, 
aniline-water and furfural-water systems was by 
measurement of the refractive index using re- 
fractive index/concentration graphs which had 
been experimentally determined beforehand. An 
Abbé refractometer was used so that only very 
small samples were needed. A _ refractometric 
method was also used for the analysis of the 
water phase in the ethyl acetate-water, ethy!] 
acetoacetate-water and isobutanol-water systems, 
while the water contents of these solvent phases 
were determined by Kark Fischer titration. It 
was found convenient to standardise the Karl 
Fischer reagent by titrating a saturated solution 
of the solvent after each experiment. The con- 
centration of ethy! formate in water was obtained 
by hydrolysis of the ester, and this method was 
also checked by cloud point titration. Formic 
acid in benzene was determined by titration with 
N/10 NaOH and the same method was used for 
the formic acid phase itself. 


Calculation of Mass Transfer Coefficients 


It was assumed that equilibrium existed at the 
interface and that the rate-controlling step in the 
saturation of either phase is the rate of transfer 
of the other component away from the interface 
into the bulk of that phase. If equilibrium exists 
at the interface then the interfacial concentration 
is constant and substantially equal to the 
saturation value. 
We then have the rate of transfer given by :** 


°*Ilam indebted to Dr. J. T. Davies for the supply of 
this compound. 
*° The symbols used are defined at the end of this paper. 
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j=k(e, —c) (1) 


where ( 
force. 


¢,, — ¢)is normally taken to be the driving 
Now the transfer rate is also given by : 
Ve (2 
A dt 
V de 
A (¢,, c) dt 


2-303 V d 
A dt 


log — (3) 


Hence log (c c) was plotted against time and 
@ 

in all cases good straight lines were obtained from 

the slope of which k could be calculated. Since 

both phases were sampled, two graphs were 

plotted for each run from which the two individual 

transfer coefficients could be obtained. Typical 


plots are shown in Fig. 2. 


@ water 
? FURFURAL 
6 
4 
2 a 
> WATER 
Et acetate 
6 
4 
2 b 
waree 
x BUTAN 
oF 
c 
™ 
Fig. 2. Typical plots of log (¢,, — ©) v. time. (4) Furfural 
~ Water at 50°C, (6) Ethyl Acetate - Water at 20, 
(c) Isobutanol — Water at 20°C. 
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Results 


The first two systems used, furfural-water and 
ethyl] acetate-water, were studied in much greater 
detail than the subsequent systems since a wider 
range of Reynolds numbers could be used before 
the interface was disrupted. On plotting the 
coeflicients against Reynolds number on logarith- 
mic co-ordinates using the Reynolds number of 
the other phase as parameter a family of smooth 
curves were obtained (Fig. 3). Occasionally very 
low coefficients were obtained and it was assumed 
that these were due to surface contamination. In 
order to check on this most experiments were 
duplicated. With the exception of points of this 
nature, the scatter was within 10°, for furfural- 
water and 20°, for ethyl acetate-water. 
Because of the large number of experimental 
results obtained with these two systems the data 
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Fig. 3. Plots of individual transfer coefficients v. Reynolds 
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plotted in Fig. 4 are smoothed values taken from 
the graphs. With all other systems the actual 
experimental data are shown. 

With a few exceptions all experiments were 
carried out with the stirrers contra-rotating. In 
order to ascertain whether the direction of stirring 
was important a number of experiments were 
conducted with the stirrers running in the same 
direction. There was no marked difference in the 
transfer coefficients apart from a tendency for the 
coefficients to be slightly greater. 


INTERPRETATION OF RESULTS 

Correlation of Transfer Coefficients 

Due to the complexity of the problem no funda- 
mental theoretical approach was possible, and an 
attempt was therefore made to correlate the data 
using dimensional analysis. Thus the mass 
transfer coefficient for phase 1 may be assumed 
to be expressed in terms of the variables for 
phase 1 and 2 as follows : 


k, =f (L, N,, Ny Pi» Par Nar D,,, (4) 


This takes into consideration the experimental 
conditions that the stirrer lengths are both equal 
to L, and that the ratios of L to the cell diameter 
and to the cell height are constant. 

Since there are 10 variables on the R.HLS. of 
(6) by the w theorem there will be (10-8) i.e. 7 
dimensionless groups. The simplest arrangement 
of these is as follows : 


(LN,*) 
Dy 


i.e. 


k,/v, =f, [(na/m), (Ap/p,), (Re,), 
(Re,), (We,) (Fr,), (Seh,)) (5) 


The values of the various dimensionless groups 
were calculated for the present systems but with 
the exception of the Reynolds groups and the 
viscosity ratio none of the groups on the right 
hand side appeared to have any relevance ; thus 
for any given system, Ap and y are constant, but 
widely different values of k/v were obtained for 


the two phases when the other groups were 
approximately equal. If Ap is not relevant then 
g can probably be neglected also since these terms 
should appear as a product with the same ex- 
ponent. This means that the second, fifth and 
seventh groups in equation (5) vanish. Again, 
the Schmidt groups varied from 100 to 14,000 
without any particular correspondence to values 
of k, indicating that this group was not effective. 
Further evidence that this was the case was ob- 
tained from experiments at temperatures up to 
50°C when only a slight increase in the transfer 
coefficients was obtained, at the same Reynolds 
numbers as at 20°C, despite the large increase in 
diffusion coefficient. 

We have thus found experimentally that 


k,/v, = fy (Re,, Rey, 


On plotting the transfer coefficients against 
Reynolds number on logarithmic co-ordinate 
paper using the second Reynolds number as 
parameter it was found that both Reynolds 
numbers were important in determining k, (cf. 
Fig. 3). From the experimental results it was 
evident that the Reynolds numbers should appear 
as a sum and not a product. A large number of 
possible relationships were tested, the dimension- 
less group k,/v, being plotted against various 
combinations of the two Reynolds numbers. The 
best correlation was obtained when plotting : 


against (Re, + Re, "2/) 


all the coefficients, with the exception of those 
for the furfural phase falling on a curve with a 
maximum deviation of about 30%, Fig. 4 (a). 
This curve can be represented, with a slight loss 
of accuracy by the equation : 


60k, /v, = 6-76 x 10° (Re, + Re, ,/n,)'** +1 (6) 


and a plot of this is given in Fig. 4 (b). 

There was a tendency for the data relating to 
any one phase to lie at a shallow angle to this 
correlation indictaing that some unknown factor 
had been omitted. The data also tended to 
diverge more from the correlation when the two 
Reynolds numbers were widely different. 

The furfural phase coefficients are plotted in 
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Correlation of individual transfer coefficients. (a) k,/v,; plotted against (Re, + Rez 72/7), 
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Fig. 4(b) and it is clear that they do not corres- 
pond at all to the correlation. No obvious reason 
for this discrepancy could be found but since the 
corresponding aqueous phase coefficients correlate 
well it cannot be due to any macro interfacial 
effect. At 30° and above the data for furfural fall 
in line with the ethyl acetate and water results. 
The coefficients obtained at higher temperatures 
i.e. 30°-50°C were also plotted in the above 
manner, Fig. 5. These data scattered more than 
at 20 
to increase with temperature. 


but there was a definite tendency for k/» 
When allowance 
is made for the decrease in v with temperature it 
can be seen that & itself remains within about 
+ 40%, of its value at 20 

Some of this divergence may be due to the 
the 
determining step were a diffusional or a chemical 


up to 50 


experimental scatter of results. If rate 
one then the coefficients would increase several 
fold over this temperature range. 

Although the the 
transfer rates, degree of turbulence and physical 


exact relation between 
properties of the two phases is undoubtedly very 
complicated, it is clear from the lack of depend- 
ence on the Schmidt group that transfer takes 
place by eddy and not by molecular diffusion 


over the range of Reynolds numbers considered. 


Effect of Interfacial Films 

Since the diffusion coefficient is normally related 
to the viscosity it was considered of interest to 
investigate the effect of changing one without the 
other. For this purpose collodion was added to 
the solvent phase in the ethyl acetate-water and 
furfural-water systems in order to increase the 
viscosity markedly without altering the diffusion 
coefficient appreciably. Unfortunately a film of 
collodion was deposited at the interface. It was 
found that although the solvent phase coefficients 
could be predicted satisfactorily from the corre- 
lation, the predicted water coefficients were far 
too low and to ascertain whether the presence of 
the rigid interfacial film of collodion altered the 
transfer rate, experiments were carried out with 
a rigid protein film deposited at the interface. 
Here there was no change in bulk phase properties 
but a very large increase in the interfacial 
viscosity. 


The mechanism of mass transfer of solutes across liquid-liquid interfaces — I 


It was found that the results both with the 
protein systems and with the collodion system 
were best explained by the assumption that rigid 
interfacial films largely prevent transfer of turbu- 
lence across the interface. Hence the second term 
in the correlation, Re, 7,/,, representing the 
contribution of the turbulence of the second phase 
to the turbulence in the first phase, is best omitted 
With 


the collodion systems the second term is small 


although a small residual effect remains. 


when correlating solvent phase transfer since 


N. <7, and its retention makes little difference. 


In the water phase, however, the second term is 
much larger than Re, and its retention thus leads 
to erroneously high predicted values. The results 
are compared with the correlation in Fig. 6. The 


Fig. 6. Effect of interfacial films on correlation. 


4 
| q 
a 
‘ 
7 \ OL ; 
10) ] 54 
20 
MOBILE 
Film 
] 
= j 
= 
ELA water! A 
$0 + o 
40 7 
RIGIO t 
10 PROTEIN + t 
1000 4000 10000 1000 4000; 10000 | 
2000 6000 20000 2000 e000 20000 


J. B. Lewis : The mechanism of mass transfer of solutes across liquid-liquid interfaces — I 


experiments with mobile adsorbed films on the 
interface indicated that there was no significant 
change in transfer rates with one molecule of 
quaternary ammonium compound per 150A?, but 
on increasing the concentration fourfold a marked 
increase was then obtained. This coincided with 
a noticeable change in the wetting angle at the 
baffle. Since the interfacial tensions were reduced 
to very low values the absence of any marked 
effect is further evidence of the non-dependence 
of transfer rates on interfacial tensions. 

In general these experiments agree with the 
hypothesis that the transfer mechanism is con- 
trolled by eddy not molecular diffusion. The effect 
of rigid films at the interface suggest that the 
%/, factor should be replaced by some inter- 
facial function possibly involving the interfacial 
Contradictory results reported by 
the effect 
agents on interfacial transfer rates may be due 


viscosity. 
various authors on of surface active 
to the marked difference in effect of mobile and 


rigid films. 


Disc USSION 


The most surprising feature of the present results 
is that the data for the transfer of water into 
furfural fail to correlate with the remainder of 
the results. All attempts to bring the furfural 
transfer coeflicients into line with, for example, 
the ethyl! acetate coeflicients, caused the data for 
The 
the pendant drop 


other systems to deviate. furfural water 


interface was examined by 
method [12] to see if saturation was accompanied 
by any marked interfacial turbulence. The effect 
appeared to be small and less than that observed 
in the other solvent water pairs used. In any case 
such turbulence, unless it gave rise to an emulsion 
in the solvent phase only, should influence the 
The 
existence of a micro emulsion on the solvent side 


transfer coefficient in the water phase also. 


cannot however be altogether ruled out. It is 
interesting to note that at 30° and above the 
furfural coefficients fall in line with the correla- 
tion. A possible explanation for this is that 
polymerisation occurs at this temperature, the 
viscous polymers formed inhibiting micro emulsi- 
fication and reducing the transfer rate to the 


normal value. 


The finding that the diffusion coefficient has no 
influence on the rate of transfer is of great 
importance, but it is not yet clear whether this 
is true with other forms of liquid-liquid contactor 
since Murpocu and Pratt have recently shown 
that the solvent-phase mass transfer coefficients 
for wetted wall columns vary approximately as 
the two-thirds power of the Schmidt group [13]. 
However, the present correlation gives the purely 
physical aspects of the transfer which in some 
respects is analogous to the collision frequency 
in gas reactions except that here the process is 
bulk, 
If slow chemical stages are involved 


due to eddy, ie. and not to molecular 
transport. 
at the interface the rate will depend on the 
energy and entropy of activation, but since the 
temperature- 
independent energy must be 
small. The standard free energies of transfer 
were calculated from the activity coefficieints at 
infinite dilution using derived from 
solubility data [14]. Table 2. It can be seen that 
there are marked differences in AG, between 
corresponding phases. Although the activation 
energies would of course be different it might be 
expected that they would follow the same trend 
as AG,. Since, however, both phases fit the 
correlation it is clear that no chemical step can 
The value of AG, for furfural is 
similar to that for the other solvents. 

The effect of entropy changes on the transfer 
rates cannot be tested directly but since, with 
one exception, all the phases fit the correlation 
then either the entropy changes are the same or 


transfer coefficients are almost 


the activation 


values 


be involved. 


else they are zero. 

It was anticipated in view of the addivity of 
the Reynolds numbers for the two phases in Fig. 4 
and equation (6) that the mass transfer coefficients 
would be lower when the stirrers were running in 
the same direction than when they were contra- 
rotating. However, experiments in which the 
stirrers rotated in the same direction indicated 
that the coefficients were, if anything, slightly 
higher than previously. The most likely explana- 
tion is that transfer to the interface is governed 
by the vertical component of the mixing length, 
which being perpendicular to the direction of 
rotation is governed by the absolute and not the 
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Table 2. Activity coefficients at infinite dilution and standurd free energies of saturation. 


(a) Aqueous Systems 


Ethyl Ethyl 


Water log fy 1-80 1-75 
Phase | 


| Acetate Formate Acetoacetate Furfural Aniline 


Ethyl 


(b) Non-aqueous System 


| Benzene | Aniline 
log fy 1-05 0-87 
AG, | 1-40 1-16 
Formic Acid | Cyclohexane 
log fy | 1-28 1-21 
AG, | 1-71 1-62 


relative values of the Reynolds numbers. * With 
contra-rotating stirrers a zone will exist where 
the mean horizontal velocity is zero but in general 
this will occur at some distance away from the 
interface in one phase or other depending upon 
the relative stirring speeds and viscosities of the 
two phases. Although this might cause a 
diminution of the turbulence in that plane the 
effect at the interface would be smal! and 
probably swamped by the eddies set up at the 
sharp edges of the baffles. The net result is 
therefore that the transfer coefficients would 
probably be almost independent of the direction 
of stirring but possibly slightly higher when the 
stirrers are rotating in the same direction. 


CONCLUSIONS 

(1) A new type of transfer cell has been designed 
and used for the determination of individual 
mass transfer coefficients for eight binary 
systems, in six of which one of the phases was 
water. 

(2) The data for all systems, with the single ex- 
ception of the transfer of water into furfural, 
have been correlated graphically by means of 

a plot of k,/v, against the group (Re, + Re, 


or 


8 


%2/,)- This can be approximated by means 
of equation (6). 


(3) Over the range of stirring speeds used the rate 
of transfer appears to be independent of mole- 
cular diffusion coefficient, and the transfer 
therefore takes place entirely by eddy mixing. 
The interfacial tension and density difference 
of the phases also appear to be without 
effect. 


(4) The effect of rigid films at the interface is to 


suppress the transfer of turbulence from one 
phase to the other so that the term in the 
correlation involving Re, drops out. On the 
other hand, the effect of mobile films at the 
interface is negligible, except at the higher 
concentration where the transfer rate is some- 


what increased. 


Most of the experimental work was carried out 
by Mr. E. W. Suarratr. The author wishes to 
thank Mr. A. S. Wurre for his interest in the 
work and Dr. H. R. C. Pratt for much useful 
advice particularly on the interpretation of the 
lata. 


NOMENCLATURE 
A area of interface in transfer cell — em*. 

¢ concentration — gm/100 cm? at time ¢. 
€ concentration at infinite time, i.e. at equilibrium. 
D diffusion coefficient — cm? /sec. 

function. 

fo activity coefficient at infinite dilution. 

g acceleration of gravity — cm/sec*. 

j transfer rate — gm/cm? sec. 

k individual mass transfer coefficient — cm/sec. 
I 
N 


length of stirrer — cm. 
stirring rate — revs. per sec. 


q 


q 
q 
Solvent log fy 1-02 1-17 0-85 0-84 0-73 0-74 5 
Phase | AG, 186 1:56 | | ove 
| 
1-76 1-82 1-32 1-68 a 
2-35 2-43 176 28 ik 
VOL, 
3 
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time — secs. 
V volume of cell compartments — cm*. 
AG, _ standard free energy of saturation — K.cal/mole. 


LN? 
Fr Froude Group -——. 


y interfacial tension — dyne/cm. Sch Schmidt Group - a 

viscosity — gm/sec. cm. D 

vy kinematic viscosity — — cm*/sec. We Weber Group 
p density gm/cm?. 


Subscripts 
L?Np 1 phase under consideration. 


Re Reynolds number 
7 2 other phase. 


Dimensionless Groups 
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The mechanism of mass transfer of solutes across liquid-liquid 
interfaces 


Part II.—The transfer of organic solutes between solvent and aqueous phases q 


J. B. Lewis 


Chemical Engineering Division, Atomic Energy Research Establishment, Harwell, Didcot, Berks. 
( Received September 1954) 


Summary—*“ Overall” transfer coefficients were determined for a number of organic ternary 
systems and the results compared with the values predicted assuming zero interfacial resistance 
and using individual transfer coefficients calculated from the correlation derived in the preceding 
paper. 
acetate-water with acetone, gave larger transfer rates than predicted, whilst others, especially 
toluene-water systems gave lower results. The faster rates were found to be associated with 
interfacial turbulence. 


Although good agreement was obtained in many cases, some systems, notably ethyl 


The low results are ascribed to slow heterogeneous reactions at the interface, and a method 
is devised for deriving the rate constants. Experiments were carried out at different temperatures 
in order to obtain the energies of activation. The probable reason for the slow step is discussed. 


Résumé— Pour un certain nombre de systémes ternaires organiques, on détermine des coefficients 
de transfert globaux qu’on compare aux valeurs prévues en supposant nulle la résistance a 
linterface et en utilisant les coefficients de transfert individuels calculés 4 partir de la formule 


établie dans l'article précédent. 


Quoique la concordance obtenue soit bonne dans de nombreux 
cas, certains systémes, notamment acétate d'éthyle 


eau, avec l'acétone, donnent des vitesses 


de transfert plus importantes que prévu tandis que d'autres, particuli¢rement toluéne-eau, 


donnent des résultats plus faibles. 


Les vitesses les plus élevées sont dues A turbulence a l’inter- 


face. On attribue les faibles résultats a des réactions hétérogénes lentes a interface et on imagine 


une méthode permettent de déduire les constantes de vitesse. 


On fait des expériences a diverses 


températures en vue d’obtenir les énergies d’activation. On discute la raison probable d'une sor 
marche lente. 
In Part I [1] the factors influencing individual is thus possible to predict the individual co- = i. 
transfer coefficients for the mutual saturation of _ efficients for solute transfer in immiscible solvent- 7 » 
binary solvent-water systems were examined. water systems where the mutual saturation - : 
Under the experimental conditions used mass technique is not applicable. a 
transfer from the interface takes place by eddy In any three-component system the overall a 
not molecular diffusion. The individual coefficients transfer rate of the undistributed component will er 
were correlated by plotting k,/v, against [Re, + depend on the two individual coefficients and on Va 
Re, n2/7,) (Fig. 4 (a) Part 1). This correlation the rate of transfer across the interface itelf. By ; a 
was represented with a slight loss of accuracy calculating the overall coefficient, assuming oe 
by the expression : infinitely fast interfacial transfer, and comparing Bean a 
” the result with the experimentally determined as) 
60 k,/v, = 6-76 x 10° [Re, + Res ng ma] ra coefficient, the value pom interfacial rate con- 4 


Since the transfer coefficients were found to be 
independent of the molecular diffusivities this 
correlation will also apply to the transfer of an 
undistributed solute in solvent-water systems. It 


stant can be determined. This will depend on 
the nature of the reaction occurring at the inter- 
face, i.e. solvation, dimerisation, etc., and hence 


information can be obtained about this step. 
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In the present investigation overall transfer 
coefficients were determined for a number of 
organic ternary systems and compared with the 
values predicted by assuming zero interfacial 
resistance. From the difference in the two co- 
efficients the interfacial transfer coefficients were 
calculated. 


EXPERIMENTAL 
Systems Used 


The apparatus employed was the same as that 
described in Part I. Six solvent-water pairs were 
used, the solvents being ethyl acetate, butyl 
acetate, tetrachloride, 
“ chlorex (bis 8 chloroethyl ether), and toluene, 
with 
benzoic acid and diethylamine as solutes. 


iso-butanol, carbon 
aniline, 
The 
distribution coefficient data for these systems 
were determined experimentally and are sum- 
marised in Table 2 and Figs. 3-9, while the 


acetone, acetic acid, ammonia, 


relevant physical properties of the solvents are 
given in Table 1. Diffusion coefficients were 
determined for all the systems but since it was 
demonstrated in Part I that these are not needed 
for the correlation the results will be reported 
elsewhere. 

The experimental procedure was similar to that 
used in the case of the mutual saturation experi- 
ments, except that solute was added initially to 
one of the phases. Both phases were sampled as 
before, but since somewhat larger quantities were 
needed for analysis samples were only taken every 
ten minutes. In a few cases, as, e.g., with dilute 
solutions of acetic and benzoic acids in toluene 
and carbon tetrachloride, very large samples were 
needed to obtain accurate analyses and inter- 
mediate samples were not therefore taken. In 
most experiments the phases were not mutually 
saturated at the start of the experiment, although 
where the solubility was appreciable, as for ethy| 
acetate, iso-butanol and butyl acetate, a few 
duplicate runs were carried out with the phases 
initially saturated. However, no significant 
difference in behaviour was observed. A few 
transfer runs were also carried out with mono- 
layers of a surface active agent adsorbed at the 
interface, the long chain quarternary ammonium 


salt C,,H,,N(CH,),Cl being used as in Part I. 


Analysis 

The acids and bases were determined by titra- 
tion with N/10 or N/40 NaOH or HCl, using a 
potentiometric method for the determination of 
the end point. (The latter was particularly valua- 
ble in the case of the solvent phase as most 
indicators are soluble in the solvents used.) 

Acetone was determined by the hydroxylamine 
method [2] using N/10 NaOH for the titration, the 
end point being determined potentiometrically in 
all but the first few runs. Aniline was determined 
by bromination in the usual way, using KBrO,- 
KBr mixture, followed by titration of the excess 
iodine with N/10 Na,S,O, after addition of KI. 


Determination of the Experimental Overall Mass 
Transfer Coefficients 

The transfer of a third component from phase 1 
to phase 2 involves transport up to the interface 
in phase 1, transfer across the interface, and 
transport away from the interface in phase 2. As 
the accumulation at the interface due to adsorp- 
tion must be negligible compared with the bulk 
concentrations the rates of these three steps will 
be equal. The concentrations on either side of 
the interface are normally assumed to be in 
equilibrium. This assumption is discussed further 
below. Since the interfacial concentrations cannot 
be determined it is customary to define an overall 
transfer coefficient based on the assumption that 
the resistance to transfer lies entirely in one 
phase [3]. If, for example, the overall coefficient 
is based on the water phase, K,,, the assumption 
is made that there is no concentration gradient in 
the solvent phase the interfacial concentration 
being equal to the bulk value c,. The interfacial 
concentration on the water side, c,*, is taken to 
be in equilibrium with c, so that for transfer from 
solvent to water we can express the rate as : 


j = K,, — ¢) (1) 


Although, in general, this assumption gives a 
purely imaginary picture of the process the 
correct transfer rates can be calculated if K,, 
is known. If the distribution coefficient H is 
constant ¢,* will be equal to c,/H. When, 
however, H varies with concentration a graphical 
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Table 1. Physical properties of phases. 


Density 


Buty! acetate 


Carbon tetrachloride 


Chiorex (dichlordiethy] ether) 
Ethy! acetate 


iso-butanol 


Toluene 


Aqueous solutions 
(i) 4% Acetic acid 


(ii) 7% Acetic acid 


(iii) 12% Acetic acid 


(iv) 4% Acetone 


(v) 4% Diethylamine 


(vi) 4% Aniline 


0-88 


1-60 
1-57 


1-22 


0-90 


0-80 
0-80 


0-86 
0-86 
0-84 


1-00 
1-00 
0-99 


Viscosi 
gm/mi. Centipoise | 


0-72 


0-97 
084 


2-38 


0-46 


The effect of adding 5% of 
various solutes to these sol- 
vents was investigated. It 
was found that the viscosities 
and densities remained sub- 
stantially unchanged. 


The densities and viscosities 
of the other aqueous solutions 
are substantially the same as 
pure water. 


method of determining ¢,*, or H, must be 
used. Thus, referring to Fig. 1 (a), the curve AIB 
represents conditions at equilibrium and the point 
O the bulk phase concentrations c, and c,. The 
assumed interfacial concentrations c,* and c, are 
represented by point B on the curve and the 
corresponding value of H is then the slope of the 
chord AB. As the bulk concentration varies 
during an experiment B moves along the curve, 
producing a family of chords such as AB, cor- 
responding to the change in H. Samples are 


taken at intervals so that the change of ¢, and c,, 


and therefore c,*, with time can be determined. 


Now at any time the transfer rate is given by :t 


V de,, 
A dt 
Vi dew 


A dt “(c,* —<.) (2) 


When H changes considerably during the course 


+ The symbols used are defined at the end of this paper. 
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Fig. 1. Relation between the individual and overall 

resistances when the interfacial transfer is immeasurably 

fast, (a) transfer from solvent to water, (b) transfer from 
water to solvent. 


of an experiment it is necessary to integrate this 
expression graphically. Thus we have : 


The right hand integral is evaluated by plotting 


1 
es — against c, and measuring the area 
(Cy, Cw) 
under the curve. Thus, mean values of K 


ow 
relating to successive values of c, can be obtained. 
In most cases H was either constant or varied 
little during the experiments and it was then more 
convenient to proceed as follows : 
If c = (¢, +¢,) is the total initial concen- 
tration (since the volumes of the two phases were 


equal), 


Vi 
Att’ —c,)/H — 
V Hit dey 


AH +1 —<, 


= — +1)—c,) (2a) 


Thus K,,, can be obtained from the slope of the 
curve obtained by plotting (c,/(H+1) — c,) 
against time on semi-logarithmic co-ordinates, 
using the appropriate values of H corresponding 
to the concentrations of the samples taken and 
H If is inde- 
pendent of concentration a linear plot will be 
obtained. 

For transfer from water to solvent an exactly 
analogous equation can be derived thus : 


At 

and K, can be obtained by graphical integration 
as before. When, as was usually the case, H 


taking a mean value of 


remained substantially constant during an 
experiment the expression was _ integrated 
directly, giving : 

K.. In (c,, — ¢o/(H +1)) (2b) 
where H is now the slope of the chord AB in 
Fig. 1 (b) and the appropriate values are used as 
before. 

Similar equations are obtained when using the 
overall coefficient based on the solvent phase, i.e. 
K,,; in this case the equilibrium concentration 
at the interface is c,*. 

The experimentally determined transfer co- 
efficients are summarised in Tables 3 to 7. 


INTERPRETATION OF RESULTS 
Theoretical 


(a) Immeasurably Fast Interfacial Transfer 


As stated earlier, the transfer of a solute from, 
e.g. solvent to water, involves three consecutive 
processes. When the transfer across the interface 
is immeasurably fast the concentrations on either 
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side of the interface will be in equilibrium, i.e. 
Cyi/Ce; = H. If the molecular species are different 
in the two phases as, for example, acetic acid 
between toluene and water where the acid is 
monomeric in water but almost entirely dimeric 
in toluene the true interfacial equilibrium differs 
from that determined experimentally since the 
former refers to the same species in both phases 
whereas the latter refers to the total concentration. 
Since however the concentrations used in the rate 
equations are the total concentrations, the inter- 
facial equilibria considered must be those referring 
to the total concentrations, i.e. H must be the 
experimentally determined value. Thus we have : 


J =k, — Cy) = he (Cos — (8) 


Now we normally define an overall transfer 
coefficient K,, equation (1), and it can be 
shown from equations (1) and (3) that when H 
is constant. 


(4a) 


When H varies with concentration a graphical 
method of analysis is used. Thus, referring to 
Fig. 1 (a), the point J on the equilibrium curve 
corresponds to the interfacial concentrations c,, 
and ¢,,. 

chord 1B 


If we now define H’ as the slope of the 


(Sa) 


then from equations (1), (8) and (5) 


1 
“i * Hi, 


(4b) 


The position of J can be obtained if k, and k, 
are known, since : 
Slope OF = — 


Cwi “we 


ky (from equation 3) 
k, 


Fig. 4 (a) and equation 6 of Part I enable k, and 
k, to be predicted for any solvent-water system, 
and the overall transfer coefficient can therefore 
be predicted whether H varies or not. 

During the course of an experiment the bulk 
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phase concentrations will vary, point O moving 
to some point O, at time ¢t. Since the volumes of 
both phases are the same the increase in con- 
centration of the water phase will equal the fall 
in concentration of the solvent. Hence line O O, 
will have a slope of — 45° and the concentrations 
of the phases during any experiment will move 
along this line. For every point O, there will be a 
new value of the interfacial concentration J, and 
O., will be parallel to OJ. Since the slope of OI 
is generally not widely different from that of O O, 
the interfacial concentrations will change much 
less than the bulk concentration, i.e. while O 
moves to O,, I only moves to J, and B to B,. 
Since H’ is given by the slope of 1B the change in 
this factor depends on the shape of the curve 
between J and B. Apart from the toluene-water 
and carbon tetrachloride-water systems the change 
in H’ was negligible in the present work. 

For transfer in the reverse direction we have : 


Referring to Fig. 1 (b), points O, J and B refer to 
concentrations in the bulk, at the interface and 
to the hypothetical overall interfacial concen- 
tration c,* and c, respectively. Thus defining 

Cui 

1 1 4 1 


(4c) 


The position at J is obtained as before and will 
again change during each experiment. 


(b) Slow Interfacial Reaction 


When the interfacial transfer rate is slow, e.g. in 
the case of a slow chemical reaction, the con- 
centrations on either side of the interface are no 
longer in equilibrium with one another. Since the 
reaction is heterogeneous it is convenient to 
express the reaction rate as an equivalent rate 
either on the water or solvent side of the interface. 

For transfer from solvent to water if the 
reaction is assumed to occur on the water side 
of the interface the concentration at the interface 
in equilibrium with the solvent will be c,,* =c,,/H. 
Then for a first order reaction the rate will be 
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equal to k,, (c,;* —c¢,;). In most cases the 
interfacial step involves reaction between the 
solute and one of the phases and since the latter 
will be in considerable excess we can assume that 
the reaction is pseudo-first order. 
exception to this is the dimerisation of carboxylic 
acids in aprotic solvents.) At the steady state 
there is no build up of solute at the interface 
and the three steps will take place at the same 
rate. 


(A_ possible 


he (Cui (6) 


= (¢,.* — as before. 


If H is constant, then 


and from equation (6) it can be shown that 


1 1 


rw 


Now the first term on the R.H.S. of this equation 
is equal to the reciprocal of the overall transfer 


coefficient, i.e. 1/K,,, giving : 
1 1 1 1 
* Hk, * k (7a) 


ow rw 


The sum of the first two terms on the R.H.S. of 
equation (7a) is equal to the reciprocal of the 
overall coefficient for the case where the reaction 
is immeasurably fast, i.e. to 1/K,..;.4;.)- Hence 


substituting for the latter in equation (7a), 


1 1 a 
k 


rw K (obs.) (ealc.) 


(8) 


If the reaction is not of the first order, it has been 
shown [4] that for an nth order reaction equation 
(8) becomes : 

(Co* — Ce) 1 1 


K ww (obs.) K (cate.) 
An analogous argument applies when the direction 
of transfer is from water to solvent. 

Instead of considering the overall water film 
coefficient, K,., the overall solvent film co- 
efficient, K,,, may be used instead. Similar 


expressions are derived and the reaction is now 


assumed to occur on the solvent side of the 


interface. For a first order reaction we have : 
_ _ 
k, K,, (obs.) K,, (cale.) 
and for a reaction of the nth order : 
rs s s 0s (obs. os (calc.) 


It is frequently found that the reaction would 
be of the first order if occurring in one phase but 
of a higher order in the other. In such cases 
therefore it is convenient to consider the overall 
coefficient, and hence the reaction, in terms of 


b 
~ 
A 
Cw Cw Cw 


Fig. 2. Relation between the individual and overall resist- 

ances in the presence of a slow interfacial reaction, (a) 

transfer from solvent to water, (b) transfer from water to 
solvent. 
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Cs 
| sO 100 200 250 300 330 400 
20 
30 
Fig. 8. Distribution of aniline between toluene and water at 20°C. 
70 
60 
40F 
t= 
Ol i. = i i 
100 iso 2 750 300 400 450 $00 $50 
Cw 
8 
? 
+f 
re) i i i i 
270 40 60 8 10 120 140 160 
Cw gn/t 
Fig. 4. Distribution of acetic acid between toluene and 
water, (x 20°C, 30°C, A BC). 
7 
20F 
10} 
Cw — 
Fig. 5. Distribution of acetone between toluene and water ‘Fig. 6. Distribution of diethylamine between toluene and 
at 20°C, water, (© 20°C, A 380°C). 


the phase which would give first order kinetics 
as this simplifies the calculation and interpreta- 
tion of results. Pratr and Murpocu for example 
[4] considered the overall coefficients for urany] 
nitrate transferring between solvents and water 
in terms of the solvent phase, as this gave first 
order kinetics, whereas basing the overall! coeffi- 
cients on the water phase gave a third order 


reaction. 


When H varies with concentration it is again 
necessary to resort to a graphical method of 
analysis. Referring to Fig. 2(a), curve AKJB 
represents the conditions at equilibrium and O 
the bulk phase concentrations. The point J 
representing the-two interfacial concentrations 
C,.; and ¢,, will not lie on the curve since the two 
interfacial concentrations are not in equilibrium. 
If we define H* as the slope of the chord JB, 
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22 " — 60) 
a bh 
40h 
or 
e 
+ > 20} 
| 
‘ 
| f A A. 
| 20 40 60 80 
7 | Fig. 9. Distribution of acetone between CCl, and watre 
| at 20°C. 
| 4h 
~ Table 2. Systems with constant distribution coefficients. 
iy 
| Solvent — Water Solute H=c,/¢, | Temp. °C 
o4F 2 | 
/ — 
| Isobutanol-Water | Acetic acid 1-18 20° 
0 a} ro Isobutanol-Water | Ammonia 0-26 20° 
Isobutanol-Water | Diethylamine | 8-88 20° 
4 Isobutanol-Water Diethylamine 5-30 80° 
; : Water | Acetone 0-81 20° 
Fig. 7. Distribution of benzoic acid between (a) ccl, and Butyl acetate- 
water, and (b) toluene and water (© 20°C, 30°C, 50°C) Water | Diethylamine 1-00 20° 
Ethyl acetate- 
c c . ° 
H* (9a) Water | Acetic acid 0-84 20 
Co” — Coi* Ethyl acetate- 
Water | Acetone 1-24 20° 
then it is found from (6) and (9a) that 
. (6) (9a) Chlorex-Water Acetone | 1-20 20° 
A Toluene-Water | Diethylamine | 0-84 380° 
c c..) — + 
" 
k, k, H*k, 
60 
50 
E 30- 
20F a 
° $0 100 iso 200 250 300 3$0 400 450 soo $s0 
Cw 
10 
LOW CONCENTRATION of © 
b+ 
=> 4) 
se 


° 20 40 60 80 100 120 140 160 180 200 220 240 
Cy 


Fig. 8. Distribution of acetic acid between CCl, and water at 20°C. 
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1 
kre K ow (oded.) (cale.) 
In order to obtain the second term on the R.H.S. 
of equation 7 (b), i.e., the predicted value of K,,, 
it is necessary to know the position of the point J 
on Fig. 2 (a) so that H* can be determined. Since 
I does not lie on the intersection of OJ and the 
equilibrium curve the position of J and therefore 
of J is uncertain and a series of approximations 
must be made. The slope of OJ is equal to 
— k,/k, as before, and if OJ is projected to meet 
the curve at L we can taken the first value of 
H*, H,* to be the slope of BL. Using H,*, the 
first value of k,, is calculated, k,.,. Now 


i.e. (7b) 


J = (Cop? — Sui) 


~ 


Since the term (c,,* — ¢,,) corresponds to JJ, 
using k,, a first approximation to (c,.* — ¢,,) 
is calculated, giving by simple geometry the first 
position of J and J. Thus a second value of H*, 
H,* is obtained from the slope of J,B and the 
procedure repeated until a constant value of k,, 
obtained. 

During an experiment point O moves along 
line O 0, due to concentration changes, as before, 
so that both H and H* change. These variations 
must be taken into account if the curvature of 
AKJB is very unsymmetrical between J and B. 

For transfer in the reverse direction an 
analogous situation applies as shown in Fig.2 (b). 
Defining H”° as follows : 


H° = (9b) 


i.e., aS the slope of the chord BJ, we can show 
that equation (8) again applies. A similar series 
of approximations may be necessary to obtain 
the correct value of H?. 
can be expressed in terms of the solvent phase 
if more convenient. 


Again the coefficients 


Application to Present Data 
of Calculated and 


(a) Comparison Observed 
Coefficients 
The experimentally determined overall water 


coefficients K,, are summarised in Tables 3 to 8 


along with the values calculated assuming im- 

measurably fast interfacial reaction. The ratios of 

K 

_owlobsd:) should equal unity if the assumption 
ow (cale.) 

of immeasurably fast reaction is correct. Values 


in the range of, say, 0-8-1-2 may well be obtained 
however because of the uncertainty in k, and 
k, due to the scatter in the correlation and to 
experimental errors, particularly the difficulty in 
avoiding interfacial contamination in the transfer 
cell. Although many results fall in this range it 
is clear that in some experiments the observed 
transfer coeflicients are much greater, and in 
others much smaller than the predicted values. 
The interfacial reaction constants k,, have been 
calculated where the observed coefficients are 
significantly lower than the calculated values. 


(b) “ Normal” Transfer Coefficients 
Good agreement is obtained between the observed 
and calculated overall coefficients with transfer 
in the following solvent-water systems : 
Tsobutanol 
(i) Acetic acid from isobutanol to water. 
(ii) Ammonia from water to isobutanol. 
Butyl Acetate 


Acetone transferring in both directions. 


Chlorea 


Acetone transferring in both directions, 


Carbon Tetrachloride 
(1) Acetone from water to CC],. 


(ii) Acetic acid from CCl, to water. 


Toluene 
(i) Acetic acid, diethylamine and acetone 
from water to toluene (there is a tendency 
however, with the last two solutes for 
low values to be obtained at the highest 
teynolds numbers, this is discussed fur- 
ther below). 


(ii) Aniline from toluene to water. 


(c) High Transfer Coefficients 
Significantly high transfer coefficients are obtained 
with transfer in the following systems : 
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394 


410 


20 
20 
20 


20 
20 


Reynolds Numbers 


W ater 


2,400 
2,400 
2,400 
2,400 
2,600 
2,600 


3,000 
3,000 
3,000 
3,000 
3,000 
3,000 


2,800 
2,800 
2,800 
2,800 


3,000 
3,000 
3,000 
3,000 


3,000 
3,000 


8,000 
3,000 


Solvent 


1,000 
1,000 
1,000 
1,000 
1,000 
1,000 


1,000 
1,000 
1,000 
1,000 
1,000 
1,000 


1,000 
1,000 
1,500 
1,500 


1,000 
1,000 
1,500 
1,500 


1,000 
1,000 


1,000 
1,000 


Conc. gm 


Initial 


Final 


Cale. 


Iila. Diethylamine transfer 


Water 


26-1 
30-0 
56-1 
25°8 


22-6 


Solvent — Water 


> Solvent (Interfacial turbulence) 


te 


31 


IIIb. Acetic acid transfer 


Water 


| 49-8 
| 52-6 
47-8 
44:8 


48-0 
38-9 
45-4 


40°5 
41-7 
35-9 
43-6 


Solvent — Water 


38-8 
31-4 
32-7 
34-7 


1:17 
1-17 
1-64 
1-64 


1-15 
1-15 
1-69 
1-69 


IIIc. Ammonia transfer 


Solvent — Water (Interfacial turbulence) 


1-4 
1-2 


Water — Solvent 


18 
1-6 


0-9 
0-9 


0-42 
0-42 


0-42 
0-42 


Koy cm/sec > 


» Solvent (Interfacial turbulence) 


108 


Obsd. 


1-17 
2-34 
2-50 
2-66 
1-83 
1-67 


1-33 
1-33 
1-33 
1-33 
2-00 
1-67 


1-33 
1-50 
1-83 
1-83 


1-67 
1-67 
2-17 
2-33 


0-25 
0-38 


0-83 
0-50 


* Concentrations of the phase from which transfer is taking place. 


®@ Values of k,, are only tabulated where they are markedly significant. 
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Run Temp. k,, cm/sec. 
+ 38) 20 35-3 
447 20 
x 458 20 39-7 
387 30 29-4 
484 30 31-6 1-31 
| 
| 4 
380 20 27-1 25-5 | 
- 448 20 35-9 34-9 
452 20 72-3 68-6 
4 459 20 34-4 32-1 
388 30 30-0 28-6 : 
485 30 31-1 30-6 
| | | 
363 
368 
364 
370 20 
362 20 
366 20 
365 20 
4 369 20 
| | 
4 
486 | 1-7 | 
423 | | 
| | | | 
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Table 4. Ethylacetate-Water : Interfacial turbulence in both directions in all systems. 


cm/sec x 108 


Temp. 
°C tial Cale. Obsd. 


Acetone transfer 
Water -> Solvent 


2,500 | | 28-1 
2,500 
4,000 155 
4,000 
5,000 19-2 
5,000 366 (195 
5,000 01 | 252 


Solvent — Water 


41-5 31-0 
33-6 
26-7 
63.2 43-3 
37-7 26-9 
41-1 20-9 


IVb. Acetic acid transfer 
Water — Solvent 


0-8 27-9 
1046 67-0 


Solvent — Water 


45-9 29-7 
90 (580 


Isobutanol In all these systems transfer was accompanied by 
(i) Diethylamine from water to isobutanol. ™«rked interfacial turbulence. This effect, which 

is related to spontaneous emulsification, was 
frequently seen when the interfacial tensions of 
Ethyl Acetate the unequilibrated three-component systems were 
determined [5]. For this purpose, the pendant 
drop method was used [6], a droplet of the solvent 
phase being suspended in the aqueous phase and 
Carbon Tetrachloride an enlarged image of the drop profile projected 
Acetone transferring from CCl, to water. on a screen. In several instances the interface 
could be seen to be in a state of violent agitation 

Toluene during transfer of the undistributed third com- 
Aniline transferring from toluene to water. ponent. Except with the two ethyl acetate-water 


(ii) Ammonia from isobutanol to water. 


Acetone and acetic acid transferring in both 
directions. 
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R. olds Numbers Conc. gm 
No. x 108 
— = = 
284 20 2,300 | 0-59 2-5 
287i 20 | 2,300 | 0-59 2-8 | 
270 | 20 8,600 1-16 3-7 
274 20 3,600 1-16 4-2 
275 20 4,600 2-01 3-7 
283 20 4,600 2-01 42 
460 20 4,600 | 201 4-0 
VOL, 
285 20 | 2,500 2,500 | 0-59 2-8 3 
7 286 20 2,500 2,500 0-59 3-2 1954 
230 20 4,000 4,000 | 1-22 27 
263 20 4,000 4,000 “7 
279 20 5,000 | 5,000 1-77 40 
885 20 5,000 5,000 1-77 42 ae 
461 20 5,000 5,000 1-77 3-7 
462 20 | 4,100 5,000 42 
438 20 | 5000 | 5000 | | 1-56 3-2 
| | | = 
a 
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Table 5. Butyl acetate-Water. 
Reynolds Numbers Conc. gm Koy cm/sec x 10° 
Run Temp. — OM / sec. 
No. *¢ Water Solvent Initial | Final Cale. Obsd. x 13 
Va. Acetone transfer 
Water —> Solvent 
307 20 2,700 3,000 0-0 | 33-8 0-73 0-5 - 
308 20 2,700 3,000 31-0 28-8 0-73 10 - 
309 20 4,500 5,000 | 320 | 210 1-68 1-9 - 
Solvent — Water 
304 20 3,000 3,000 42-6 } 28-8 0-76 0-8 - 
306 20 3,000 3,000 42-6 25°8 0-76 0-8 - 
389 20 5,000 5,000 31-7 | 21-0 1-73 12 - 
411 20 5,000 5,000 387 (256 1-78 1-5 - 
Vb. Diethylamine transfer 
Water —> Solvent 
428 20 4,000 5,000 33-6 28-7 | 1-73 0-8 15 
| 
Solvent + Water 
427 20 5,000 5,000 | 27-9 23-4 1-89 10 2-1 
Table 6. Chlorex-Water 
| Reynolds Numbers | Conc. gm/t* K,. cm/sec x 108 
No. °C Water Solvent | Initial | Final Cale. Obsd. om sec 
Acetone transfer 
| Water -> Chlorex 
492 20 «| 3,000 61-6 43-6 1-54 2-32 
494 20 2,700 3,000 | 480 1-54 1-77 
495 20 | 3,600 | 4,000 | 527 | 834 2-41 2-65 
Chlorex — Water 
491 20 | 8,000 3,000 | 462 37-8 1-62 1-82 
493 20 | 3,000 3,000 47-6 36-7 1-62 1:58 
496 20 | 4,000 4,000 42-8 | 32-3 2-55 1-82 
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Tabie 7. Carbon tetrachloride-Water. 


Reynolds Numbers Cone. gm K,,,. cm/sec x 108 


W ater Solvent Initial Final Cale. Observed 


Vila. Acetone transfer 
Water Solvent 


3,600 4,000 41-9 38-0 0-66 
3,600 4,000 36-9 34:1 0-66 
3,600 4,000 410) 39-5 0-66 
3,600 4,000 40-5 38-6 0-66 


Solvent — Water (interfacial turbulence obsd.) 


4,000 
4,000 31-2 


VIIb. Benzoic acid transfer 
Water — Solvent 


1-6 
1-7 


0-95 
Solvent — Water 


1-95 1-85 
1-03 1-87 
3-59 3-35 
20 1-00 
3-69 3-37 


Ge 


Vile. Acetic acid transfer 
Water —> Solvent 


119-5 118°8 0-038 
103-4 103-0 0-038 
121-6 120-2 0-038 
120-5 120-0 0-032 


Pax 


Solvent — Water 


28-6 
25-0 


35-2 


+ Interface probably contaminated. 


3 
Run Temp. . Kap Cm 
No. c x 108 
445 20 0-63 
+4 20 0-52 
468 20 0-35 
ere 
443 20 4,000 24-0 0-85 1-67 VOL, 
*. 446 20 4,000 25-4 0-85 1-00 3 
: 1954 
454 20 4,000 4,000 18 2-02 0-92 1-8 
= 477 20 4,000 4,000 10 1-77 0-62 1-0 tal. 
e 456 30 4,000 4,000 1-9 1-64 0-70 1-2 a 
479 30 4,000 +,000 10 1-44 0-43 0-6 
455 20 4,000 4,000 1-33 13 
469 20 4,000 4.000 0-67 10 
476+ 20 4,000 4,000 0-17 0-2 
457 30 4,000 4,000 3-04 gt 
480 30 4,000 4,000) 195 
| 
| 
| 
467 20 3,200 4,000 0-03 
465+ 20 3,300 4,000 0-02 0-04 a 
478 20 3,300 4,000 0-05 ad 
481 30 3,300 $,000 0-05 
| 
464 20 4,000 $000) O14 O-15 
466 20 4,000 40-7 0-14 0-20 
482 30 4,000 4,000 30-6 0-138 0-17 
i 
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Table 8. Toluene-Water 
Reynolds Numbers | Conc. gm /t* Koy cm/sec x 10° 

Run Temp. —-— Om /sec. 
No. °C | Water Solvent Initial Final Cale. |  Obsd. x 108 

— 

| VIIla. Acetic acid transfer 

| Water — Solvent 
371 | 20 2,800 3,000 53-2 52-8 0-083 0-050 _ 
375 20 2,800 3,000 53°38 53-0 0-037 0-033 - 
475 20 2,500 3,000 1145 (1128 0-078 0-098 
397 30 2,900 3,000 7-5 47-2 0-025 0-042 - 
418 1) 2,700 3,000 | TS 708 | 0-033 0-075 | - 
376 | 20 4,700 5,000 53-4 53-0 0-084 0-050 0-12 
379 l 20 4,700 5,000 548 54-2 0084 | 0067 | 0-88 
414 | 20 4,700 5,000 | 615 61-3 0-084 0-033 0-06 
435° 20 4,700 5,000 (40.1 0-084 0-017 0-02 
436(a)t 20 4,700 5,000 | 425 | 423 0-084 0-017 0-02 
436(b)t 20 | 4,700 5,000 42-4 | 42-2 | 0-084 0-033 0-06 
474 20 4,200 5,000 1235 | 1217 #+| ©20 0-67 mt 
398 30 4,800 5,000 | 54-2 54-0 0-062 0-025 0-04 
417 50 4,800 5,000 | 65-5 | 0-077 0-10 
487 50 4,600 | 5,000 131-0 1289 0-83 - 

Solvent — Water (Interfacial turbulence observed) 
373 20 | 3,000 3,000 459 | 333 | 0-176 | 0-183 - 
374 20 | 3,000 3,000 43-7 | 82-5 0-176 | 0-130 0-5 
396 30 3,000 3,000 43-2 33-0 O14 | 0-15 - 
416 50 | 8,000 3,000 436 | 288 0-108 | 0-167 - 
877 20 5,000 5,000 41-6 30-6 0-418 | 0-012 0-02 
378 20 5,000 5,000 450 | 834 | O418 | 0083 —§ O10 
412 20 5,000 5,000 46-5 28-5 O-418 | 0-167 0-28 
413 | 20 5,000 5,000 51-7 28-4 0-418 0-217 0-45 
434(a)f | 20 | «5,000 5,000 42-0 30-7 0-418 0-167 | 0-28 
434(b)t | 20 5,000 5,000 28°3 21-4 | 0-367 0-13 0-20 
399 30 5,000 5,000 44-7 29-2 0-343 | = 0-238 0-72 
400 30 5,000 5,000 46-1 32-6 0-343 0-183 0-39 
415 | 5O 5,000 5,000 43-6 26°5 0-23 0-20 1-36 
488 50 5,000 5,000 (569 37-6 023 | 0-25 
VIIIb. Diethylamine transfer 
| Water — Solvent 
347 | 20 1,600 2,000 32-3 302 O81 0-28 - 
313 20 1,600 2.000 39-5 37-9 | 0-31 0-18 
392 | 20 1,600 2,000 | 327 30-2 0-21 0-35 “ 
296 20 2,400 3,000 | 52-7 46-9 0-56 0-43 - 
297 20 2,400 3,000 | 403 35-9 056 | 050 - 
383 30 | 2,500 3,000 38-0 33-9 0-58 | 0-52 
489 30 | 2,500 3,000 35-6 31-4 0-58 | 0-57 - 
310 20 8,900 5,000 | 39-5 34-3 1-20 | 0-57 10 
390 20 | 8,900 5,000 32-6 27-8 1-20 1-08 71 
391 20 8,900 5,000 | 325 | 28-7 120 | 1-8 
Solvent — Water 

312 20 | 2,000 2,000 40-8 36-4 | 0-85 0-27 18 
848 20 2,000 2,000 35-5 31-0 0-35 0-27 18 
204 20 | 8,000 3,000 44-0 36-1 0-63 0-47 1-9 
295 20 | $000 | 8,000 33-8 28-2 0-63 0-37 0-9 
882 30 | 8,000 | 8,000 319 +%| 2-7 0-59 0-42 15 
490 30 8,000 | 3,000 34-2 31-2 0-59 0-52 44 
298 20 | 5,000 | 5,000 40-4 31-4 1-49 0-42 0-6 
$11 20 | 5,000 | 5,000 40-7 82-3 1-49 0-50 0-8 


+ Monolayers on interface. 
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Reynolds Numbers 


W ater Solvent 


3,000 


353 20 2.700 3.000 
346 20 4,500 5.000 
20 4,500 5.000 

4 5,000 


Solvent 


3,000 3.000 
352 20 3,000 3.000 
351 20 5,000 5.000 
5,000 5,000 


508 20 5,000 5,000 
499 5O 5,000 5,000 
5O2 5O 5,000 | 5.000 
507 SO 5,000 5,000 


509 5O 5,000 5,000 


431 20 5,000 5,000 
433 20 5,000 5,000 
500 5O 5,000 5,000 


501 5000 | 5,000 
5,000 5,000 
510 50 5,000 5,000 


5 
g 


3,000 
3,000 


20 
20 38,000 | 3,000 


| 3,000 3,000 


Ville. Acetone transfer 


>» Water (Interfacial 


VIIId. Benzoic acid transfer 


Water —> Solvent (Interfacial turbulence observed) 


a 

q 
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Table 8 cont. Toluene-Water. q 

q 


Conc. gm 


Initial Final 


Water — Solvent 


42-6 0-61 
40-3 36-7 0-61 0-60 
33-7 32-6 144 0-22 0-3 


35-2 27°8 144 1-17 6-2 
29-8 1-44 


turbulence observed) 


48-9 0-69 


31-6 29-9 0-69 0-67 
37-9 28-4 1-55 0-98 2 
29-7 24:55 1-55 0-380 1-7 


Water — Solvent 


1-79 1-46 2:38 1-17 
0-62 0-45 2-32 1-08 20 
3-51 2-65 1-52 95 2°5 
2-91 2-15 1-50 13 
1-22 O75 1-25 “38 

2-04 | 1-51 18 55 


Solvent — Water 


= 
zx 
= 
=: 
~ 
~ 
te 


2 1+ 
2-04 1-04 2-62 1-17 2 
1-64 1-38 1-59 2-62 
1-68 1-53 1-59 1-48 40 
7-37 6-85 2-97 2-33 41 


eo 
t 
te 
x 


Aniline transfer 


1.16 
1-16 


17.6 
20-8 


Solvent — Water 


1-00 6-6 


19-0 18-7 1-17 
21-6 21:3 1-17 0-83 2-8 


systems, where transfer of acetone and acetic acid 


in either direction caused turbulence, this effect 
if obtained at all was confined to transfer across 
the interface in one direction only. The direction 
of transfer causing turbulence is recorded in 
Tables 3 to 8, for those systems where this effect 
was found. 

The increased turbulence and larger area arising 
from this effect will increase the observed transfer 
rate by an unpredictable amount so that in the 
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absence of a slow interfacial reaction the ratio of 


K (obed.) 
however, a slow reaction does occur then the net 


will be greater than unity. If, 


result will be complicated and ratios equal to or 


less than unity might be obtained. Thus, in a 
system where this effect is found a value of the 


obs./cale. ratio near unity indicates the probable 


existence of a moderately slow interfacial reaction. 
When a ratio significantly lower than unity is 
obtained this is still greater than the value that 


a 
| | | 
Run Temp. - Kup q 
No. Cale. Obsd. x 108 
| 
| VOL. 
32 | 20 
2 5,000 5,000 
1954 
| 0 
| | 
| 
4 432 20 2,900 | | 13.2 1.67 
426 20 | 2,900 | 1-67 
425 
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would be obtained if allowance could be made 
for the additional turbulence and increased area. 
Consequently the value of k,, calculated will be 
too great in all systems where interfacial turbu- 
lence arises. The effect on the calculated activa- 
tion energy cannot be predicted. 


(d) Low Transfer Coefficients 

Significantly low transfer coefficients are obtained 
in a number of instances and for these the reaction 
rate constants were calculated as described above. 
Where this could be done at more than one tem- 
perature an approximate value of the activation 
energy E was obtained. In a few instances the 
rates are low at 20°C but “ normal” at higher 
temperatures so that only the 20° value of the 
reaction constant could be obtained. It was 
thus not possible to calculate the activation 
energies but either these were very small or else 
the low coefficients at 20° were not significant. A 
few transfer coefficients were also obtained that 
were ‘ normal * at Reynolds numbers of 2,000 and 
8,000, but somewhat low at 5,000. 
that here the reaction constants are of the same 
order as the turbulent transfer coefficients at 
Reynolds numbers of 5,000 but sufficiently greater 
than the turbulent coefficients at lower Reynolds 
numbers not to have a marked effect on the 


It is clear 


transfer rate. 

The transfer of acetic acid from water to 
isobutanol and from toluene to water gives rise 
to marked interfacial turbulence. In these cases 
therefore high transfer coefficients should be 
obtained and their absence indicates a slow 
reaction. With isobutanol the coefficients are 
approximately normal” so that no reaction 
constant can be calculated. Low coefficients are 
obtained with the toluene system but the calcu- 
lated reaction constant will be too large and the 
activation energy very approximate, 

The systems where slow reactions probably 
occur are: 


Isobutanol 
(i) Diethylamine from isobutanol to water 
at 20°C, 
Since obs./cale. ratios greater than 
unity are obtained at 30° either the 
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(ii) 


activation energy is small or else the low 
values at 20° are not significant. 


Acetic acid from water to isobutanol. 
Ratios slightly greater than unity are 
actually obtained but since an appreci- 
able interfacial turbulence was observed 
these values indicate the probable pre- 
sence of a slow interfacial reaction. The 
rate constant cannot be obtained but it 
is probably of the same order as the 
turbulent transfer coefficients. 


Butyl Acetate 


Diethylamine transfer in either direction. 

Since experiments were only carried out at 
20° no value of the activation energy was 
obtained. 


Carbon Tetrachloride 


(1) 


(ui) 


Benzoic acid in either direction at 20°C, 
For transfer from solvent to water at 
30° ratios greater than unity were ob- 
tained thus either E is small or the 20° 
result is of no significance. In the other 
direction the position is not clear as slow 
transfer also occurs at 30° with a similar 
reaction constant to that at 20°. This 
may well be due to experimental in- 
accuracies, 
Acetic acid from water to CCl, at 20°C, 
A ratio greater than unity is obtained 
at 30° so that again either E is small 
or the 20° results are of no significance. 


Toluene 


(i) 


(ii) 


Acetic acid from toluene to water. 

Interfacial turbulence was present in 
this system so that the observed overall 
coefficient and the calculated reaction 
constants are both too large. The appar- 
ent energy of activation is about 4 
K.cals./mole but due to the turbulence 
no assessment of the accuracy of this 
value can be made. 


Diethylamine from toluene to water. 
E is somewhere between 5 and 7 
K.cals./mole. 
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(iii) Benzoic acid from water to toluene. 
E is about 10 K.cals. /mole. 


(iv) Benzoic acid from toluene to water. 
E is about 13 K.cals./ mole. 


In addition to these results low values were also 
obtained with between 
toluene and water when the interface was coated 
with a monolayer of C,,H,,N(CH,),CIl. For trans- 
fer from toluene to water the rate was approxi- 
mately the same as in the absence of a film but 
in the other direction a marked retardation in 
transfer was obtained as compared with transfer 
across the clean surface. 


transfer of acetic acid 


Discussion 

Of 82 systems studied 15 gave reasonable agree- 
ment between the calculated and observed 
coefficients, while 8 transfer rates were faster and 
9 slower than predicted. All the fast transfer 
rates were obtained in systems exhibiting marked 
interfacial turbulence and since the 
“ normal ” and only two of the slow rates exhibited 
any interfacial effect this phenomenon is con- 
sidered to be the cause of the high rates. Since 
“normal ”’ rates were obtained, in some cases at 
least, for every solvent-water pair investigated, 
with the exception of ethyl acetate-water which 
was, however, actually used in obtaining the 
correlation, the results were considered to be a 
justification of the use of the correlation in pre- 
dicting third-component transfer in the transfer 
cell. The diminution in the rate of transfer of 
acetic acid from water to toluene caused by the 
monolayer was unexpected since a similar film 
had little effect om rates of mutual saturation of 
ethyl acetate-water and furfural-water. It may 
be however that such monolayers are very specific 
in their action. Recent work on the retardation 
of evaporation from water surfaces covered with 
unimolecular films tends to support this view [7]. 

Except for transfer between toluene and water 
at the highest Reynolds numbers, the acetone 
transfer coefficients were either in agreement with 
the predicted values or else high due to inter- 
facial turbulence. The only effect likely to occur 
with acetone transfer would be and 
entropy changes arising from the different type 


none of 


energy 


of solvation in the two phases. This should be 
most marked with toluene-water and CCl,-water 
where the solvents and water are most dissimilar. 
It is interesting to note that the distribution 
coefficient of acetone between toluene and water 
has a marked temperature coefficient, unlike dis- 
tribution between the ethyl acetate-water and 
butyl acetate-water phases. The slight retardation 
in transfer between toluene and water is thus 
probably due to such a solvation effect. Since the 
effect is small it would not be noticeable at the 
The 
barrier would of course be higher in one direction 
than the other and if this is in the direction toluene 
to water, the retardation in transferring from 
the 
analogous case of transfer from carbon tetra- 


lowest stirring speeds. interfacial energy 


water to toluene may be negligible. In 


chloride to water a marked interfacial turbulence 
is obtained which will mask any small interfacial 
effects. 

Acetic and benzoic acids are known to dimerise 
in aprotic solvents such as toluene and carbon 
tetrachloride whereas in water, alcohols, ketones, 
etc. the monomeric molecules are stabilised by 
hydrogen bonding to the solvent molecules. The 
transfer of these acids between the toluene-water 
and CCl,-water phases will thus involve a dimeri- 
the 
absent in the other systems investigated. It is 


sation reaction at interface which will be 
significant therefore that slow transfer rates have 
been obtained (with the exception of acetic acid 
from CCl, to water) for the transfer of acetic and 
benzoic acids between toluene and water and 
CCl, and water. Transfer of acetic acid between 
the other solvents was not slow, except possibly 
from water to isobutanol where turbulence may 
have masked a small effect. 

The reaction constants have been calculated 
assuming first order kinetics. This may well be 
incorrect but in the absence of other evidence this 
is probably the best way of expressing the results. 
The values of k,, and E obtained are only 
approximate as the aim of the present work was 
to examine the general behaviour of a large 
number of systems rather than obtain very 
accurate data on a limited number. (More detailed 
investigations are now being carried out on a few 


systems of particular interest.) No data are 
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available as to the free energy and heats of 
dissociation of acetic and benzoic acids in toluene 
but in benzene the values are 4-3 and 8-2 K.cals. /- 
mole respectively for acetic acid and 5-2 and 8-4 
K.cals./mole for benzoic acid [8]. Similar values 
apply for acetic acid in carbon tetrachloride. Since 
however the interfacial reaction occurs in the 
presence of water the nature of the reation will 
be different and these values of no relevance. 
Transfer from solvent to water must involve a 
reaction between dimeric acid molecules and water 
whilst transfer in the other direction involves a 
reaction between monomeric acid molecules 
hydrogen bonded to water molecules. Warp [9] 
has calculated the increase in free energy of acetic 
acid transferring from water to toluene at 25°C 
from distribution data. This is about 3 K.cals./- 
mole but no value of the activation free energy 
can be computed ; allowance has also to be made 
for the energy of adsorption. Since benzoic acid 
is much less soluble in water than acetic acid the 
degree of solvation may be widely different and 
this may alter the reaction rate appreciably as 
compared with acetic acid. Further, with in- 
creasing concentration of acetic acid in the toluene 
phase the solubility of water will increase. This 
will then reduce the extent of dimerisation in this 
phase and hence accelerate the rate of reaction, 
as was actually found. Further information as to 
the nature of the rate-controlling step is clearly 
desired. It is interesting that the turbulent 
transfer technique is able to demonstrate activa- 
tion energy barriers when the technique of 
diffusion across a plane interface fails [9]. 

Diethylamine was found to transfer slowly 
between butyl acetate and water and from 
toluene to water. A small effect was also found in 
transferring from isobutanol to water. The effect 
here is almost certainly associated with the pre- 
sence of the amine group which presumably 
forms hydroxyl compounds in water. Diethyla- 
mine being a stronger base than ammonia should 
thus have a more marked retardation of transfer. 
This appears to be the case as the observed 
ammonia transfer coefficients agree with the 
calculated value for transfer from isobutanol to 
water whereas the corresponding diethylamine 
values are somewhat low. 


CONCLUSIONS 


(1) Overall transfer coefficients were determined 
for various solutes transferring between six 
different solvent-water pairs, 16 systems, i.e. 
32 cases being investigated in all. The results 


were compared with values calculated using 
the correlation derived in Part I. In 15 cases 
the two sets of data agreed well, in 8 the 
observed rates were faster and in 9 slower 
than predicted. 


(2) The high coefficients were found to occur only 
in systems where marked interfacial turbu- 
lence arose and it was concluded that this 
phenomenon was the reason for the high 
values. 


(3) The low coefficients were associated with the 
transfer of those solutes which might reason- 
able be expected to be involved in reaction at 
the interface. Approximate values of the 

first-order reaction constant k,, were calcu- 

lated, and in a number of instances very 
approximate values of the activation energy 
of the reaction also. 


(4) It is concluded that the correlation of Part I 
is quite satisfactory for predicting individual 
turbulent transfer coefficients in the transfer 
cell used. In the absence of any interfacial 
turbulence and slow chemical reaction, the 
overall transfer coefficient can then be calcu- 
lated. Conversely if interfacial turbulence is 
absent, the chemical rate constant can be 
obtained if the observed and calculated over- 
all transfer coefficients are known. 


(5) It is thus demonstrated that the behaviour of 
two different systems should not be compared 
even in the same equipment unless it is known 
that interfacial turbulence is absent and that 
no slow interfacial reaction is present. 


The bulk of the experimental work was carried 
out by Mr. E. W. Suarratr. Thanks are due to 
Mr. A. S. Wurre for his continued interest in this 
work, and to Dr. H. R. C. Pratt for many helpful 


discussions. 
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NOMENCLATURE 


A area of interface in transfer cell — em* L_length of stirrer — cm. 
¢ concentration gm /100 cm*, at time N,, Ng stirring rates in phase 1 and phase 2 
€9 total initial concentration. respectively — revs. per sec. 
c* concentration in equilibrium with other t time — secs. 
phase. V volume of cell compartments — em*. 
E activation energy — K.cals /mole. L?N, 
H_ distribution or partition coefficient —c¢,/c,, Ke, Reynolds number in phase 
H’,H, H*,H° defined by equations (5a), (5b), (9a) and I sn 
(9b). Re, Reynolds number in phase 2 = —$ 
j transfer rate — gm/cm* sec. "2 
K,, overall" mass transfer coefficient based viscosity gm/sec em. 
on solvent phase cm/sec. vy kinematic viscosity = — em*/sec. 
corresponding overall coefficient based on p density gm/em*, 
water phase — cm/sec. Subscripts 
k individual mass transfer coefficient 1 phase under consideration. 
em /sec. 2 other corresponding phase. 
k,, velocity constant for first order interfacial s solvent phase. 
reaction based on solvent phase — cm /sec. si value in solvent phase at interface. 
corresponding velocity constant based on water phase. 
water phase — cm/sec. wi value in water phase at interface. 
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A rapid method of measuring absorption rates and its application to 
CO, absorption into partially carbonated ammonia liquor 
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Summary 


An experimental absorption system designed to simplify the investigation of pro- 


cesses in which the absorption rate is limited by liquor reaction is described. 


Data from an investigation into the mechanism of absorption of carbon dioxide into partially 


carbonated ammonia liquor is utilised to illustrate the operation of this absorption system and 


is then employed for determining the kinetic constant for the reaction controlling the absorption 


rate in this process. 
by other techniques. 


Résumé 


The values obtained agree satisfactorily with those previously measured 


L’auteur décrit une méthode expérimentale d’absorption congue pour simplifier 


l'étude des processus dans lesquels la vitesse d’absorption est limitée par la réaction en phase 


liquide. 


Ii utilise les résultats d'un travail sur le mécanisme d’absorption du CO, dans une solution 


ammoniacale partiellement carbonatée, pour illustrer le mécanisme de cette méthode d’absorption 


et calculer ensuite la constante cinétique de la réaction contrélant la vitesse d’absorption. Les 


valeurs obtenues sont en bon accord avec celles qui avaient été mesurées auparavant par d'autres 


procédés. 


INTRODUCTION 


Any thorough experimental investigation into a 
particular case of absorption accompanied by 
chemical reaction in a packed tower usually 
involves a considerable number of measurements 
of absorption rates under differing conditions of 
gas composition, liquor flow rate, liquor composi- 
The 


liquor and gas analyses which must be made 


tion and temperature. large number of 
accurately in order to obtain absorption rates 
and driving forces, together with the necessity of 
possessing comprehensive vapour pressure data 
for the liquor, are liable to deter the intending 
investigator whose time is limited. The apparatus 
the task of 


obyv lates the 


necessity of using vapour pressure data for the 


to be described greatly simplifies 


measuring absorption rates, and 
interpretation of the absorption mechanism. This 


apparatus has been used for determining the 
mechanism of the absorption of carbon dioxide 
into partially carbonated ammonia liquor, and 
utilised in 
Data 


from these measurements are also used to obtain 


results of this investigation will be 


order to illustrate its method of operation. 


the temperature variation of the kinetic constant 
for the second order liquor reaction between free 
ammonia and free carbon dioxide. 


THe EXPERIMENTAL METHOD 


A flowsheet of the apparatus is given in Fig. 1. The 
absorption tower consists of a lin. diameter 
vertical glass tube containing 76 ceramic discs 
threaded on a vertical wire running down the 
axis of the tube. Each dise is set with its face 
at right angles to the faces of the discs imme- 
diately above and below it. The gas stream 
passes up the tube and surrounds the dises which 
are wetted by the descending liquor stream. The 
tube walls are not wetted. The general arrange- 
ment of the absorption tower is similar to that 
described in Ref. [1], and is particularly suited 
to the measurements of absorption rates where 
the liquid film is controlling. The design shown 
has a gas/disc interfacial area of 355 sq. cm. 
The gas flow up the tower is maintained by a 
positive displacement rotary gas booster which 
circulates gas through the tower at a rate of about 


4 M*/hr. This corresponds to a gas velocity of 
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about 8 ft/sec up the tower. The gas leaving the 
top of the tower passes through a reservoir vessel 
fitted with a by-pass and purge connections before 
re-entering the booster. The circulating gas is 
virtually at atmospheric pressure, but is isolated 
from the atmosphere by the small lute filled with 
transformer oil shown teed in to the descending 
gas pipe. 

When the absorption tower was in use for 
absorbing carbon dioxide in partially carbonated 
ammonia liquor, pure carbon dioxide was fed into 
the circulating gas stream through the connection 
to the descending pipe at a rate that was sufficient 
to replace the carbon dioxide absorbed in the 
liquor passing down the tower. The feed rate was 
manually adjusted so as to maintain the circu- 
lating gas pressure at the lute equal to atmospheric, 
the carbon dioxide absorption rate being then 
directly indicated on the capillary flowmeter. 

As the process of carbon dioxide absorption is 
liquid-film limited, the rate of absorption is not 
increased by increasing the gas flow rate up the 
absorption tower. With the design of tower shown, 
the amount of carbon dioxide absorbed from the 
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Fig. 1. Arrangement of experimental apparatus. 


gas stream in a single pass up the tower was only 
a small fraction of the total amount passing up 
the tower due to the high gas circulation rate and 
the low absorption rate; the partial pressure of 
carbon dioxide in the gas stream decreasing by 
less than 5°, between the bottom and top of the 
tower. Variation in the liquor composition and 
temperature between top and bottom of the 
tower was also small. The absorption tower is 
thus virtually an element of packing surface of 
known area in which gas of a constant composition 
is contacted with liquor of a known composition 
and where the resultant absorption rate is directly 
measured. 

In order to promote absorption it is, of course, 
necessary for the partial pressure of the absorbing 
gas in the gas stream (p,) to be higher than its 
partial pressure in equilibrium with the liquor 
(p,). The difference in partial pressure (p, — p,) 

Ap is known as the driving force promoting 
absorption, and must be evaluated in order to 
analyse the results of absorption rate measure- 
ments. This quantity is usually obtained from 
the results of a gas analysis, which gives p,, and 
the use of vapour pressure data for the liquor in 
conjunction with measurements of the liquor 
temperature and composition which gives p,. 
This procedure is unnecessary with the apparatus 
described as the absorption train and gasmeter 
attached on the outlet side of the lute enable 
the driving force (Ap) to be read directly. Their 
use is best described by considering the circulating 
gas content of the system, firstly, when the gas 
is in equilibrium with the liquor passing down 
the absorption tower, and secondly, when carbon 
dioxide is admitted into the circulating system 
and absorption is taking place. 

In the first case, the circulating gas contains 
partial pressures of carbon dioxide, ammonia and 
water vapour in equilibrium with the liquor, 
together with sufficient air to bring the total 
pressure up to one atmosphere. In the second 
case, the circulating gas contains partial pressures 
of ammonia and water vapour in equilibrium 
with the liquor as before, a partial pressure of 
carbon dioxide exceeding the previous value by 
the driving force (Ap) and a partial pressure of 
air Jess than the previous value by the amount 
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Ap. Provided that the circulating gas and liquor 
temperatures are maintained the same by suitable 
adjustments to the heaters fitted to the absorption 
tower, gas pipes, gas reservoir and circulating 
pump, then the volume of air displaced from the 
circulating stream through the lute and into the 
gas meter is (Ap)(V) where V is the volume of 
the gas circulating system. As a certain amount 
of ammonia, carbon dioxide and water vapour 
are also displaced through the lute, a small acid 
filled ammonia scrubber and a Sofnolite packed 
carbon dioxide absorber are interposed between 
the lute and the gas meter to ensure that the 
meter only records the displaced air. Knowing 
the volume of the gas circulating system and the 
temperature of the meter and the circulating gas, 
the constant of proportionality between meter 
reading and partial pressure driving force is 
readily obtained. 

In order to measure the rate of absorption of 
carbon dioxide into a liquid of known composition 
and temperature flowing at a known rate down 
the absorption tower and under a known partial 
pressure driving force, it was thus merely necessary 
to set the liquor flow rate by the rotameter, start 
the gas booster and adjust the gas and liquor 
temperatures by means of the heaters, and then 
by admitting carbon dioxide into the system 
displace sufficient air through the gas meter to 
give the required driving force. The capillary 
flowmeter reading when the gas meter pointer 
was stationary on the correct driving force and 
the lute had ceased blowing was then equal to 


the carbon dioxide absorption rate. 


In order to check values of the driving force 
obtained in the above manner it was possible to 
isolate the gas reservoir and analyse its contents 
when the apparatus was running. In this con- 
nection it may be noted that it was essential to 
maintain the whole of the apparatus in contact 
with the circulating gas at a temperature slightly 
above the liquor temperature in order to prevent 
condensation on the walls. The gas reservoir and 
gas analysis train were placed in a gas heated air 
oven, the absorption tower and descending gas 
pipe were wound with heating elements and the 
gas booster was surrounded by steam coils. All 
measurements were made with circulating gas 
compositions and temperature for which the 
formation of ammonium carbamate was im- 
possible. 

ANALYSIS OF 

Tue MECHANISM OF ABSORPTION 
Theorectical investigations into the process of 
absorption with chemical reaction [2] have dis- 
tinguished certain limiting mechanisms to which 
actual systems approach in behaviour. The 
recognition of the particular limiting mechanism 
which best describes the characteristics of the 
absorption system under study is most easily 
effected by setting down the behaviour of these 
limiting cases in response to variations in experi- 
mental constraints. The most readily variable 
constraints are: liquor flow rate (L), partial 
pressure driving force (Ap) and concentration of 
the uncombined liquor reactant (B). The following 
table gives the criteria which characterise the three 
limiting mechanisms distinguished in [2]. 


Table 1. List of absorption criteria. 


Mechanism 


Influence on absorption rate 


Description Conditions* Of L Of Ap of B 


* For list of symbols see end of art ick 
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Of the above three mechanisms, the first is very 
similar to the simple process of absorption without 
chemical reaction, inasmuch as a negligible amount 
of reaction of the dissolving gas takes place in the 
liquid film, and consequently the rate of transfer 
is solely determined by diffusion. If the reaction 
is very slow indeed then the absorption rate will 
become limited by the rate of reaction in the 
liquor held up in the tower packing. In the limit 
this fourth mechanism results in an absorption 
rate proportional to the product of the reaction 
rate and L°*. The last term expresses the fact 
that the hold up in a fully wetted packed tower 
is roughly proportional to the one third power of 
the liquor flow rate. 

In general, an actual case of absorption with 
chemical reaction will not follow any one of the 
fixed limiting mechanisms exactly. Its behaviour 
will rather tend towards one or other of the 
limiting cases as the experimental constraints are 
varied. Thus, the second of the tabulated limiting 
mechanisms is favoured by a low driving force 
and a high reactant concentration, whereas the 
third is rendered likely by the opposite conditions. 
A theoretical analysis based on the classic film 
theory of cases intermediate between the limiting 
mechanism quoted in Table 1 is given in [2]. 
Analyses based on the unsteady state theory of 
absorption give rather different expressions for 
absorption rates from the classic theory. All 
theories are, however, agreed on the general 
features of the criteria listed in Table 1, and 
agree exactly on the expression for absorption 
rates for the second limiting mechanism, but 
there is some uncertainty as to the value of the 
index of the liquor rate effect for the first and 
third mechanisms. The indices of 0-5 to 0-9 which 
are given have been taken as covering the range 
of values found experimentally. 


Tue CARBONATION OF AmMMONIA LiqQuoR 


Laboratory investigations into the absorption of 
carbon dioxide into ammonia liquor in a packed 
tower have previously been reported by Van 
KREVELEN 4 al. [3] and recently in a wetted-wall 
tower by Sitver [4]. Both these investigations 
were concerned with dilute (less than 2-2 N) 
virtually uncarbonated ammonia liquor such as is 
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used in gas works practice. The measurements 
reported here extend the range of liquor com- 
positions up to 11N ammonia liquor, whilst 
dealing systematically with the effects of driving 
force, of liquor flow rate, of liquor temperature 
and of degree of reaction, on the carbon dioxide 
absorption rate. 

The effect of partial pressure difference driving 
force on the absorption rate for three different 
liquor compositions is shown in Fig. 2. The values 
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Fig.2. Effect of CO, driving force on CO, absorption rate. 


plotted are the direct experimental readings, the 
absorption rate being the capillary flowmeter 
readings, and the driving force being the reading 
of the specially calibrated gasmeter. The series 
of readings for one of the liquor strengths is given 
in Table 2. In making this series the absorption 
tower was run continuously with a feed liquor of 
constant composition, and the driving force 
increased at intervals by 100 m.atm. steps. Due 
to this procedure both R, the mean carbonation 
ratio of the liquor in the tower (the mol ratio of 
CO, to NH,), and the mean liquor temperature 
increased slightly between the first measurement 
at 106 m.atm. driving force and the last at 927 
m.atm. Apart from this small effect it is, however, 
quite clear from Fig. 2 that the carbon dioxide 
absorption rate is proportional to the driving 
force Ap. It is also evident that the liquor com- 
position has a pronounced influence on the 
constant of proportionality between absorption 
rate and driving force. Reference to the table of 
criteria characterising the three limiting mech- 
anisms indicates that the only mechanism likely 
to satisfy the above conditions is the second. 
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e This conclusion gains further su rt when the 
3 x ppo 
rae 3/12 is effect of liquor flow rate on absorption rate is 
studied. The results of tests made with four 
= | different liquor compositions showing the effect 
lele 3 >) . of liquor flow rate are plotted in Fig. 8. In order 
| | to reduce the measured absorption rates to a 
common basis and eliminate the influence of 
g = | liquor composition, the ordinate is taken as the 
a | giz lSisié ratio of the absorption rate measured at a given 
e | s liquor flow rate to the absorption rate measured 
g S | using the same liquor and driving force at a 
| standard flow rate, which was taken to be 
3 a | & 3 10 litres/hr. All the points obtained fall on a 
. ? smooth curve which is divisible into two parts, a 
ele steeply rising curve in which the absorption rate 
sis increases rapidly with increase in liquor rate 
4 R followed by an almost flat line where the absorp- 
| i2 tion rate increases only slightly with increase in 
§ r ie cicigig liquor rate. The sharp fall off in absorption rate 
& below about 7 litres/hr is ascribed to incomplete 
| wetting of the packing at these low liquor flow 
2 3 rates, the gas/liquor interfacial area decreasing 
£ "i> ~ 1 {[gi* as the liquor rate decreases. The slow increase in 
> 
| 
= — S = —) 
i | | / 
a. | 
tis = gi = 039 
PIP IP |S O52 
= | Ek | | 6967 O41 
1 
2 o 4 8 12 6 
3) eials LIQUOR FLOW RATE be 
; : Fig. 3. Effect of liquor flow rate on CO, absorption rate. 
& ; | Z| 274 > absorption rate with liquor flow rate when once 
| 28 | the packing is fully wetted is further confirmation 
| ; & 2 = a | 
that the controlling mechanism is the second. 
listed in Table 1. The final criterion for the 
Sie mechanism was the non-existence of any influence 
rf of liquor flow rate on absorption rate. The small 
a 
B | ? | = < 3 influence measured is probably solely due to the 
| greater interfacial area at high liquor rates con- 
. 
3 g§ g§ | E =) sequent upon the greater thickness of the liquor 
ae |e 
= | | layer on the discs. 
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Having shown that the mechanism controlling 
the absorption rate over the range of liquor 
compositions, driving forces and flow 
investigated agrees with the 
limiting type, it is then possible to determine the 
kinetics of the liquor reaction from measurements 
of the effect of liquor composition and temperature 
on the absorption rate. Reference to Table | 
shows that the order n of the liquor reaction with 
of the 
reactant (ammonia in this case) is double the 


rates 


closely second 


respect to concentration uncombined 
index of the variation of absorption rate with 
reactant concentration. Thus, if the absorption 
rate divided by the driving force is proportional 
to the square root of the total ammonia strength 
of the liquor for a given temperature and degree 
of reaction, this indicates that the reaction is 
first order with respect to ammonia. Two such 
plots of experimental values for virtually un- 
carbonated ammonia liquor (carbonation mol 
ratio circa 0-05) at temperatures of 20-7°C and 
60-5°C are given in Fig. 4. 
lines are obtained up to a liquor strength of 11 N 


Satisfactory straight 


which was the highest used. The rate controlling 
absorption reaction is thus first order with respect 
to ammonia, as was also shown by VAN KReEVELEN 
et al. and Sitver, for weak liquors [3, 4 

Now for the absorption of a component 4 from 
the gas stream into a liquor in which it reacts 


CARBONATION RATIOS 0-05 TO 0 O7 
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Fig. 4. Effect of strength of liquor on CO, absorption rate. 
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with component B, under conditions such that 
the absorption mechanism is of the second type, 
and the liquor reaction has a velocity given by 


A = kAB, it has been 


the kinetic equation — 


shown [2] that the absorption rate per unit area 
of interface is : 


Absorption rate N 


4 =P,Hy VRBD (1) 


where H, = Henry’s law constant for solution 

of A 

D = diffusion coefficient of A in liquor 

k = second order kinetic constant 

p, = partial pressure of carbon dioxide in 
gas 


B = concentration of reactant in liquor 
In the present case this gives : 


Absorption rate of carbon 
dioxide per unit 
interface N, 


H, VkD(NHs) . p, (2) 


where (NH) is the concentration of free dissolved 
ammonia in the liquor. 

Equation (2) may be re-arranged so as to enable 
the kinetic constant to be calculated from the 
absorption rate measurement, thus : 


N, 


l 
DH (2a) 


The last term in the above equation may be 
obtained from the slope of the lines in Fig. 4 
divided by the interfacial area of the disc tower, 
which was 355 cm®. In consistent units values for 
the two lines are: 

9-8 10° mol /cm? sec atm (mol /em®)! at 20-7°C 


and 


20.5 < 10° mol /em?® sec atm (mol /em*)! at 60-5°C 


Taking values of 


D = 1-82 x 10° em?*/sec at 20-7°C, 
D = 4-53 « 10° cm?/sec at 60-5°C, 
and H, = 3-85 x 10° mol/cm® atm at 20-7°C, 
H, = 1-55 x 10° mol/cm® atm at 60-5°C, 
then 3-56 10° cm*/mol sec at 20-7°C, 
k = 38-6 x 10° cm*/mol sec at 60-5°C, 


From these two values of the kinetic constant its 
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variation with temperature is found to be given by 
the equation : 


log. = 14-28 — 


Values from this equation are found to agree well 
with those directly measured for the liquor 
reaction CO, + NH, NH,COOH by Faurno vr 
[5] and more recently by Pinsent, Pearson and 
Roventon [6]. 

For partially carbonated ammonia liquors it is 
to be expected that Equation (2), in a modified 
form, will similarly correlate absorption rates. 
The ammonia concentration will be the uncom- 
bined ammonia content of the liquor, not the 
total content, and the carbon dioxide partial 
the 
p,) as the partially carbonated liquor exerts 


pressure driving force will be difference 
(P, 
a partial pressure of carbon dioxide of magnitude 


p, The modified form of equation (2) is : 


N, = H, VkD(NHs) (p, — (3) 


As the reacted carbon dioxide combines in more 
than one form with ammonia and an equilibrium 
between the various ions is set up in the liquor, 
the relationship between the free ammonia con- 
tent, the total ammonia content and the carbona- 
tion ratio of the liquor is complex. In these cir- 


cumstances the uncombined ammonia concen- 
tration is best calculated from vapour pressure 
data for the liquor using Henry’s law; whence 


equation (3) becomes : 


N H, kD Hypa (P, 


P,) (Ba) 


where H, = Henry’s law constant for NH, 
P, = partial pressure of ammonia over 


liquor. 


For dilute liquors, where the solubilities and the 
little affected by the 
amount of carbon dioxide dissolved in the liquor, 


diffusion coefficient are 


it follows from equation (3a) that : 
N. 
Ap 
The applicability of this law is demonstrated 
in Fig. 5 where the effect of carbonation ratio on 
the absorption rate at constant total ammonia con- 
centration is plotted according to equation (4). 
At higher total ammonia concentrations (> 2 


oc Vv Pn (4) 


mol NH, /litre) appreciable deviations occur from 
the simple law of equation (8), and equation (8a) 
must be used, making appropriate allowances for 
the lower diffusion coefficient and lower solubilities 
in liquors containing much carbon dioxide. 

With the above qualifications it appears that 
the mechanism of absorption corresponds closely 
with the limiting case defined by equation (8a) for 
a range of liquor compositions from 0-25 to 11 N 
on total NH, and up to carbonation ratios of 0-6 
for the range of temperatures of from 20°-60°C 
and driving forces of up to 0-8 atm and for the 
whole range of liquor flow rates possible in the 
experimental absorption tower. In this con- 
nection it is interesting to note that even with 
the very different hydraulic conditions obtaining 
in a wetted wall tower, Sitver found a similar 
mechanism to operate. It is, therefore, likely that 
for similar liquor and gas compositions and 
temperatures, the controlling mechanism will be 
the same in any type of packed absorption tower. 
rate measurements made in the 
experimental tower are, therefore, directly 
applicable to the design of large scale packed 
towers in these circumstances. 
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A = Concentration of absorbed com- 


ponent in liquor 


B = Concentration of reactant in 


liquor 


C, = Concentration of total CO, in 


partially carbonated ammonia 
liquor 


C,, = Concentration of total NH, in 


ll 


partially carbonated ammonia 
liquor 

Diffusion coefficient of absorbed 
component in liquor 

Henry's law constant for solu- 
tion of A 

Henry's law constant for solu- 
tion of CO, 

Henry's law constant for solu- 
tion of NH, 

Kinetic constant for the absorp- 
tion reaction (2nd order) 

Liquid film matter transfer co- 
efficient 

Liquor flow rate 


Sreruens, E. J. and Morais, G. A.; 
Van Kreveven, D. W. and Horrizer, P. J.; Rec. Trav. chim. Pays-Bas 1948 67 563. 
Van KreEvVELEN, D. W. and van Hooren, C. J.; Rec. Trav. chim. Pays-Bas 1948 67 587. 


NOTATION 
mol/em*® 

Ny 
mol/cm* 

N, 
mol/em4 = 
mol /em3 

Pr 
cm? /sec 

Ap = 

mol /em? atm R 
mol/em® atm t 
= 
mol /cm* atm Z 
em®/mol see (NHg) 
cm/sec R 

cm? /sec 
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Order of the liquor reaction with 
respect to the liquor reactant 
Absorption rate per unit of in:ter- 
face of absorbed component 
Absorption rate per unit of inter- 
face of CO, 

Partial pressure of CO, in gas 
stream 

Partial pressure of CO, in equili- 
brium with liquor 


= Partial pressure of ammonia in 


equilibrium with liquor 


= Partial pressure driving force 


Ratio C,/C,, the carbonation 
ratio of the liquor 

Time 

Absolute Temperature 

Number of molecules of B react- 
ing per molecule of A 
Concentration of free dissolved 
in the liquor 

Refers to gas volumes at 20 
pressure. 
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Equilibrium calculations 
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Summary—aA simplification is suggested of the Poncnon and Savarir difference-point 
construction method for distillation at constant molal overflow by using a reflux-concentration 
diagram. This method is particularly applicable to dilute solutions and for graphical tray-to-tray 
analysis in azeotropic, extractive, and multicomponent separations. The difference-point 
construction is common to most unit operation equilibrium calculations. 


Résumé— Dans le cas de la distillation 4 flux moléculaire constant, l'auteur simplifie la méthode 
graphique de PoncHon et SAvanrir en se servant d'un diagramme reflux-concentration. Cette 
méthode est applicable en particulier aux solutions diluées et pour l'analyse graphique “ plateau 
a plateau " dans les séparations azéotropiques, extractives et des systémes complexes. 

Cette construction est commune aux calculs d’équilibre de beaucoup d’opérations fondamentales. 


EQUILIBRIUM CALCULATIONS 
Numerous graphical procedures have been pro- 
posed for calculating the number of equilibrium 
stages in countercurrent mass transfer operations. 
In the field of distillation one of the classical 
methods is that of Poncuon [10] and Savarir [13], 
using an enthalpy concentration diagram. This 
method allows for varying molal overflow and 
so takes into account heats of mixing, varying 
molal latent heats, etc. 

The method can only be used, however, in the 
limited number of cases in which the requisite 
enthalpy data are available. The well-known 
McCabe-Thiele method overcomes this difficulty 
by making the assumption, which is approxi- 
mately valid in a number of systems, of constant 
molal overflow. The object of this paper is to 
suggest a graphical construction, using a reflux- 
concentration diagram similar to that of the 
Ponchon-Savarit method, but incorporating the 
assumption of constant molal overflow. 

The use of a reflux concentration diagram 
involving the difference-point construction is 
useful in the stripping of very dilute solutions and 
for tray-to-tray analysis in ternary and multi- 
component distillation. Interest in this method 
of calculation also arises in the teaching of 
chemical engineering since the difference point 
construction is common to several unit operations, 


i.e. adsorption, absorption, liquid extraction, 
leaching, ion exchange, etc. 


GRAPHICAL REPRESENTATION 


In Fig. 1 two horizontal lines are drawn at 
any convenient distance apart. The lower line 
is used to represent the liquid (X) and the 
upper line to represent the vapour (Y). If 
the liquid is at the boiling point, mark off on 
the (X) line the compositions of the bottoms, 
feed and distillate, namely X,, X, and Xp. 
Draw a vertical line from X, and mark off the 


reflux ratio fA 


Join MX, and extend to the 


1 
Og 
KV bo 
Xr ro 
“4 
v4 
ar 
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line through X, at Z. The difference point Z 
is located on the line through X, by the reflux 
ratio below the feed tray. Starting from the 
top of the column the tie-line AD represents 
the top plate. The line MD is now drawn, 
and cuts the Y-line at C, representing the vapour 
leaving the second plate down. The line CF 
represents this plate, while E represents the 
vapour rising to it from the third plate, and so 
on. A similar construction, using lines drawn 
through Z, gives the composition of vapour 
and liquid leaving plates below the feed. 
When the feed is a vapour, join MY, and 
extrapolate to cut the X, line at Z’. The 
actual location of the feed point on the vertical 
line from Y, to X, is determined by g, the total 
heat needed to convert one mole of feed into 
saturated vapour divided by the molal latent 
heat. If the feed is a vapour, i.e. g = 0 and the 
feed point is Y,. If the feed is a liquid at the 
boiling point, i.e. g = 1, the feed point is at X,. 
If the feed is partially vapourised g wil! be between 
zero and 1 with the feed point being a fraction q¢ 
of the distance from Y, to X,. If the feed is a 
cold liquid g will be greater than 1 with the feed 


point below X, on the vertical line Y, to X,, 
while if it is a superheated vapour the feed- 
point will be above Y,. 

It can be shown by simple geometrical con- 
siderations that the construction gives the same 


results as the McCabe-Thiele diagram, with 
straight operating lines. It may be noted that the 
X- and Y-lines need not be horizontal as long as 
they are parallel. The condition for constant molal 
overflow is that the vapour and liquid lines on the 


enthalpy-composition chart should be straight 


and parallel. 


BENZENE-TOLUENE SEPARATION 


Fig. 2 shows the graphical construction for 
determining the number of equilibrium stages 
when the feed composition at the boiling point 
is 830-0 mole °%, benzene, the overhead composition 
is 95-0 mole %, benzene, and the bottoms com- 
position 5-0 mole % benzene. For a reflux ratio 
of 8:1 the number of theoretical trays is 12-7. 
The McCabe-Thiele construction also gives 12-7 
theoretical trays. 


Os 
Fig. 2. 


StripPinGc oF Ditute So.utTions 


For the stripping of alcohol or phenol from water 
the mole fraction of component remaining in the 
water is often very low indeed, i.e. 0-00001-0-0001 
mole fraction. This low concentration makes it 
almost impossible to enlarge the McCabe Thiele 
graphical plot sufficiently to determine the number 
of equilibrium stages in the column. To overcome 
this difficulty various workers have suggested 
the use of a log-log plot to facilitate the calcula- 
tions. This log-log construction suitably spreads 
the equilibrium curve but makes the location of 
the resulting curved operating line somewhat 
difficult. 

A simpler and probably more accurate graphical 
procedure than the log-log plot is suggested for 
determining the number of theoretical or equili- 
brium stages in the column. The procedure, 
involving the difference point construction, is 
best illustrated by the following example : 

It is desired to strip an alcohol feed of approxi- 
mately 0-08 mole fraction of ethyl alcohol in 
water to a final concentration of 0-00004 mole 
fraction ethyl alchohol. The number of theoretical 
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Fig. 3. 


plates are required for A reflux ratios of 2.0 


and 3-0 respectively. 

Firstly a vapour (¥Y) and liquid (X) line is 
drawn on log-rectangular paper as shown on 
Fig. 3. Equilibrium relationships are shown 
plotted by dotted lines. For most systems at 
very low concentrations the equilibrium ratio is 
constant and consequently as seen on Fig. 38 
the equilibrium tie lines are parallel. The vapour 
(Y) and liquid (X) lines can be placed any 
convenient distance apart. 

Starting from X,, the still composition, 
the equilibrium composition Y, is evaluated. 
Using equation (1) the solid line on Fig. 3 relates 
the liquid composition (4) to the vapour com- 
position Y,. After the first one or two trays 
at the base of the column the solid lines, or 
operating lines, become parallel thus leading 
to a simplification of the graphical construction. 


From the general operating line equation 
L’X,.., = V’Y,, + WX, (1) 
=VWV’Y+(L' — 


if Xt > 1 — -.. then equation (1) reduces to 


X, L 
Ym = = constant (2) 


Thus for reasonable reflux ratios the operating 
lines become parallel when X,,,, > X,. This 
conditio® is usually fulfilled after one or two 


trays if X,,, itself is very small. 


From Fig. 3 the number of theoretical stages 
is 6-0 for yy 2-0 and 8-5 for p= 3-0. 


REctTIFYING SECTION 


The rectification of components of low relative 
volatility often requires a diagram to be plotted 
on a large scale. This is well illustrated by 
Ropinson and GILuiLanp ([12] p. 210) in the 
purification of alcohol. The scale on Fig. 4 
can be conveniently enlarged and when the reflux 
ratio is high the following construction is some- 
times useful. 


0-9 
Fig. 4. 


The vapour (Y) and liquid (X) scales can be 
placed any distance apart with the dashed lines 
AD, CF, EH, etc., representing the tie equili- 
brium lines. The operating lines are located by 
the following relationships 


AC GE_ RR 
BD FD HF R+1 
The number of equilibrium stages is determined 


as before by zig-zagging between the equilibrium 
and operating lines. 


(3) 


TERNARY MIXTURES 


Most ternary graphical methods for equilibrium 
calculations are based on the binary key com- 
ponent McCase-THIeLe construction or use the 
triangular graph difference point construction, 
i.e. [1], Warre [18]. CHamBers [4] 
uses McCaspe-TuHIe_e constructions for ternary 
tray-to-tray analysis but the overall plots are 
somewhat confusing. 

The difference point construction used on 
Fig. 1 lends itself to tray-to-tray analysis in 
ternary mixtures. A single graph is constructed 
for each component. 

Consider the example of Ropinson and 
GILLILAND ([12] p. 219) for the system benzene- 
toluene-xylene. The feed composition is 60-0 
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mole % benzene, 80-0 mole % toluene, 10 
mole % xylene, the overhead product 99-5 mole %, 
benzene, 0-5 mole % toluene and the bottoms 
composition 0-5 mole percent benzene, 74-4 
mole °%, toluene and 25-1 % xylene. The tray-to- 
tray analysis construction is shown on Fig. 5. 
On Fig. 5c for the xylene component it is to be 


XYLENE 


os 

mol FRACTION BENZENE 
a 


Fig. 5a. 


‘ . / 


mol FRACTION XYLENE 
c 


Fig. 5c. 


noted that the change over from the stripping 
to the rectifying operating line occurs at the 
composition of xylene on the feed tray and not 
at the composition of xylene in the entering feed. 
Such constructions as shown on Fig. 5 help to 
give a clear picture how each component distri- 
Fig. 5b. butes itself throughout the column. 


os 
mol FRACTION TOLUENE 
b 
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Fig. 6 gives the tray-to-tray analysis for 
example 1 quoted by CHamuers [4]. A feed 
of 50-0 mole °% methyl! alcohol, 20 mole %, ethyl 
alcohol, 30-0 mole % water enters the column 
to give an overhead distillate of 92-0 mole % 


— 
05 ro 


mol FRACTION METHYL ALCOHOL ——= 
a 


Fig. 6a. 


Os 


mol FRACTION ETHYL ALCOHOL 
b 
Fig. 6b. 


methyl! alcohol, 8-0 mole °%, water and a bottoms 
composition of 8-0 mole °%, methyl! alcohol, 32-0 
mole % ethyl alcohol, 60-0 mole % water. 
CuamBeEeRS has expressed the ternary vapour 
liquid equilibria on Y-X diagrams for methyl 
alcohol-water, ethyl alcohol-water, with the 
third component as parameter. This representa- 
tion of the equilibrium data for azeotropic and 
extractive distillation systems permits the rapid 
use of equilibrium data in the graphical con- 
structions outlined in this paper. 


T 


“OS 


mol FRACTION WATER 


Fig. 6c. 


From Fig. 6 it will be seen that the feed enters 
tray 7. The composition entering the tray is a 
mixture of the feed and the composition leaving 
tray 8 and corresponds to 62.0 mole % methyl 
alcohol, 21-5 mole % ethyl aleohol and 16-5 mole % 
water. The composition leaving the feed tray 
is in approximate agreement with the composi- 
tion leaving tray 7. In the location of the feed 
tray in this example, it is important that there 
should be the correct percentage of water on 
the feed tray since the amount of water influences 
the equilibrium relationships and thus the 
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S. R. M. Exus;: 
obtaining of the correct overhead and bottoms 
compositions. 

The representation of tray-to-tray analysis 
as shown on Figs. 5 and 6 can also be applied to 
azeotropic and multicomponent distillation. In 
the case of azeotropic distillation the reflux 
ratio in the rectifying section is often high, 
due to the presence of the entrainer and con- 
sequently the use of the construction shown on 
Fig. 4 for above the feed tray avoids the difficulty 
of the rectifying section difference point being off 
the graph paper. 


Gas ABSORPTION 


SHERWOOD [15] uses the McCabe-Thiele construc- 
tion, whereas Brown [3] and Epmisrex [5] use 
the absorption factor graphical! plot for calculating 
the number of equilibrium stages. 

The difference point construction has not been 
widely applied to gas absorption equilibrium 
calculations although Ross [13] has used this 
construction on triangular graph paper. However, 
this method is not particularly suitable for low 
concentrations of solute in the gas stream. 

The same graphical constructions as here des- 
cribed for distillation can also be used for gas ab- 
sorption equilibrium calculations. Gas and liquid 
composition lines are drawn parallel, the composi- 
tions being expressed as lbs. of solute per |b. of non- 


| 
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solute. Fig. 7 shows the plot for the absorption 
of ammonia from air into water. The gas enters 
at a concentration of 0-1775 lbs of ammonia/Ib. 
of air and the water leaves at a concentration of 
0-125 lbs. of ammonia/per |b. of water. The 
ratio of water entering to gas leaving is 1-35. 


Lieutp EXTRACTION 


The Poncnon and Savarrir difference point 
construction has been extensively used in liquid- 
liquid extraction calculations. Nasu and HunTER 
[8] have used triangular graphs whereas THreLe 
[16], Raspat and Loncmarton [11], 
[17], [7] and other 
workers 

EicGin [6] and Brown [2] have used similar 


MaLoney and Scuvupert 


have used rectangular co-ordinates. 


graphical constructions for leaching. 


Ion Excuance Units 


Perirrer [9] has proposed a difference point 
construction method for evaluating the number 
ol equilibrium stages in a continuous counter 


current ion-exchange unit. 


NOTATION 


mole fraction in the liquid phase 

mole fraction in the vapour phase. 
moles of liquid reflux from the condenser. 
moles of liquid reflux from tray to-tray above the 
feed. 
moles of liquid reflux below the feed. 
moles of vapour above the feed. 
moles of vapour below the feed. 
moles of distillate. 
moles of bottoms. 
difference point above the feed. 
difference point below the feed. 

A, B,C, D, E, F, G, points of intersection of the operating 
line cn the vapour and liquid lines. 

Xp, X;, Xy, subseripts refer to distillate, feed, and 
bottoms. 

q = total heat to convert one mole of feed into saturated 

vapour divided by the molal latent heat. 
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Cuiayton’s The Theory of Emulsions and Their 
Technical Treatment. Fifth edition by C. G. Sumner. 
Pp. viii + 669. J. & A. Churchill Ltd., London, 1954. 
72s. net. 


For all who are interested in the theory and the practical 
applications of emulsions, the appearance of a new edition 
of CiraytTon’s handbook must necessarily be a matter of 
interest. Earlier editions have been remarkable for the 
thoroughness with which they have covered the literature 
of the subject. ‘“‘ Clayton ” has never been a book for the 
uninitiated, since it accepts most authors at their own 
valuation. Nevertheless, it has been, and remains, the 
leading book in English, and probably in any language, 
as a work of reference, since very few relevant publications 
are not mentioned in its pages. 

The present edition has been prepared by Dr. C. G. 
Sumner, for many years a collaborator of Dr. CLayTon 
and himself a notable contributor to the literature of the 
subject. The book has been brought well up to date and 
now includes chapters on emulsion polymerization and 


Book review 


CORRECTION 


A convenient method for the evaluation of vapour-liquid equilibria 
of binary mixtures 


by W. Jost and H. Rock 


On page 18, right column, two addition symbols have been inserted in formula (9) 
instead of multiplication symbols (cf. the corresponding formulae on page 20). 


paint emulsions which provide adequate introductions to 
extensive new fields of literature. 


The most serious fault of the book, persisting in the 
new edition, lies in the mode of presentation of the 
theoretical material. This is distributed over the first 
300 pages and constitutes in effect a complete text-book 
of surface chemistry, with exhaustive lists of references to 
literature often of doubtful relevance to the main theme 
of the book. It is difficult to conceive that the average 
reader, seeking a succinct account of current emulsion 
theory, could remain unconfused by this mass of material, 
most of which is not subjected to critical evaluation by 
the author. 


All this, of course, does not reduce the value of the 
treatise as a work of reference, and we must be grateful 
to Dr. Sumner for the meticulous character of his revision. 
It is to be hoped that, in a future edition, he will discard 
the superfluous and attempt a critical approach to the 
theory of emulsions. 

D. F. CHEESMAN. 
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On page 23, this volume, the correct reading for Z at 35°C (middle of the page) is : 
0 


18708; 1-228; 1:1150; 
0-9886 ; 0-9271; 09172; 
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10448 ; 
0-9102, 


; 00-9700; 00-0522; 


5 


The figures given under the heading Z are actually those for 2. Comparison shows that it 
0 0 


4 


is not reasonable to carry the approximation beyond Z. 


4 
The missing number in J at 25°C is 0-9333, in J at 55°C 0-9500 and in J at 55°C 0-9423 
0 0 0 
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